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ABSTRACT 



We have investigated S-wave bound states composed of three identical bosons in- 
teracting via regulated delta function potentials in non-relativistic quantum mechan- 
ics. For low-energy systems, these short-range potentials serve as an approximation 
to the underlying physics, leading to an effective field theory. 

A method for perturbatively expanding the three-body bound-state equation in 
inverse powers of the cutoff is developed. This allows us to extract some analyti- 
cal results concerning the behavior of the system. Further results are obtained by 
solving the leading order equations numerically to 11 or 12 digits of accuracy. The 
limit-cycle behavior of the required three-body contact interaction is computed, and 
the cutoff-independence of bound-state energies is shown. By studying the relation- 
ship between the two- and three-body binding energies, we obtain a high accuracy 
numerical calculation of Efimov's universal function. 

Equations for the first order corrections, necessary for the study of cutoff de- 
pendence, are derived. However, a numerical solution of these equations is not at- 
tempted. 
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CHAPTER 1 



INTRODUCTION 



Approximations abound in physics. Classical mechanics is merely an approxima- 
tion that works well at large distances and small speeds. If distances become too 
short, we must revert to quantum mechanics. If the speeds become too fast, we enter 
the realm of relativity. If the distances are short and the speeds are fast, then the 
approximations supplied by quantum mechanics and relativity are no longer valid and 
we must combine the two to get quantum field theory. Even quantum field theory, 
which is the basis for the Standard Model and nearly all of particle physics theory, is 
likely only an approximation. It may break down at even higher energies and need 
to be replaced by something else like string theory. 

Just because an approximation is not correct everywhere does not mean it is useless 
and should be discarded. In fact, the opposite is true. Since it is correct somewhere, 
it should be embraced. Approximations provide a way of isolating and controlling 
our ignorance. Identify only the necessary details and throw the garbage away! The 
hard part, of course, is identifying what is necessary so we do not throw away the 
baby with the bath water. 

The heart of this thesis is approximations. We analyze the low-energy quantum 
mechanical three-body problem using delta function potentials. For sufficiently low 
energies, the particles cannot discern the details of the interaction. Therefore, any 
suitably-adjusted, short-range potential should provide an accurate approximation to 
the true underlying potential. All knowledge of the true potential is characterized 
by a few parameters which can be fit to experimental data. This allows us to study 
general behavior that is valid for any low energy system. 

I shall first provide several examples of approximations in areas ranging from 
classical electrostatics to effective field theories. This will hopefully illustrate the 
basic concepts used throughout the rest of our work. Readers familiar with these 
ideas may wish to skip ahead to Section 11.41 where we discuss in more detail the 
problem considered in this dissertation. 
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1.1 Multipole Expansion 



The electrostatic potential created at some point r by a charge distribution is given 
by the formula 



V(r) 



1 



4vre 



d 3 f" 



P(f) 



i.r 



where p is the charge density of the distribution. Suppose that all of the charge is 
contained within some sphere of radius I (r' < I), and that the point at which we 
calculate the potential is very far away (I <C r). In this case, we can expand the 
denominator in the integral 
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1 



cos(fl) 




cos 2 (fl) 
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where 9 is the angle between the vectors r and f . This leads to approximations for 
the potential that result from truncating the multipole expansion pQ: 



V{f) 



4vre 



d 3 f' p(f f ) H — 2 / d 3 r~*' r' cos(8) p(f') + 
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The first term in the expansion is called the monopole term. It represents the potential 
that would be created if all of the charge was concentrated at one point. For large 
distances, it makes sense that the distribution would look like a point charge, and 
the first term reflects this. However, if the total amount of charge is zero, then 
the monopole term is also zero. Yet the charges must still create some potential. 
This potential is approximated by the second piece called the dipole term. It is the 
potential that would be created by a dipole at f' = 0. The full potential is built 
from all of these terms. It behaves somewhat like a monopole, and somewhat like a 
dipole, and somewhat like a quadrupole, etc. 

Notice that as we move farther away from the charges, all of the terms decrease 
in strength, but some decrease faster than others. For a non-zero total charge, the 
dominant term is the monopole term at sufficiently large distance, so it is called the 
leading order term. The dipole term is then called the next-to-leading order term, or 
equivalently the first order correction to the leading order term. If the total charge 
is zero, then the dipole term is the leading order term with the quadrupole moment 
providing the first order correction. 

Even if we stay at a fixed radius, the higher terms in the expansion still contribute 
less and less. Consider the dipole term. The integration involves r' which we know 
to be less than or equal to I. We expect that the integral would be roughly equal to 
I times some charge. The whole term then looks like l/r 2 . This is smaller than the 
leading term, which behaves like l/r, by a factor of l/r <C 1. We say that this dipole 
term is of order 0(l/r) compared to the leading term. 

The small quantity l/r acts as an expansion parameter for the potential. Each 
additional term is smaller and smaller. If we desire some accuracy e in our calculation, 
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we only need to include a finite number of terms from the expansion. Let N be the 
integer such that (l/r) N < e. Then any term in the expansion of 0((l/r) n ) for n > N 
may be dropped. For example, if N — 2, we keep only the monopole and dipole 
terms. 

The expansion hinges upon the fact that I and r are widely separated length 
scales. If we considered distances where r ~ /, then each term of the expansion would 
be about as large as all the others. They all contribute equally, and the expansion 
breaks down, reflecting the fact that we are now considering distances short enough to 
discern the details of the charge distribution. Thus, there is a limit imposed on how 
small r may be. Exceed this limit, and the approximation is worthless. Even without 
the r > / restriction, we find that the individual terms in the multipole expansion 
diverge as r — > even if the true potential never diverges. This provides additional 
proof that the multipole approximation is no good at very short distances. 

1.2 Quantum Mechanics 

The Schrodinger equation in position space is 

iP(r) = E^(r). (1.4) 

We will choose units so that H — 1. Let us consider the case of a spherically symmetric 
potential and look at low-energy S-wave scattering from this potential. We will not 
cover scattering theory in great detail, but rather try to treat the subject quite simply. 

If the potential is zero, then solutions to the equation for E > are easily found. 
There are two linearly independent solutions which we take to be incoming and out- 
going spherical waves. Any S-wave solution is written as a linear combination of these 
two solutions: 

pifcr _ — ikr 

rf}(r)=A + B . (1.5) 

r r 

The energy for this wavefunction is E — k 2 /2m. For a non-zero short-range potential, 
this solution will still be valid at large distances where we can neglect the interac- 
tion. So we can view scattering as a spherical wave approaching the potential from 
r = oo, interacting with the potential, then leaving the potential and returning to 
infinity. Since the probability associated with the incoming and outgoing waves must 
be conserved, the constants A and B can only differ by a phase, A = —Be 2l5o( - k \ The 
function So(k) is called the phase shift, and the k dependence is shown to illustrate 
the fact that the phase can depend upon the energy of the wave. 

Analytic properties of the phase shift So(k) for very low momentum can be used 
to show that the quantity k cot(<5 (A:)) has a well-defined expansion about the point 
k — 0. This expansion, known as the effective range expansion, is written as 

k cot(5 (fc)) = -- + ~r e k 2 + ■ ■ ■ . (1.6) 
a 2 
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For a more detailed derivation of this expansion, see Ref. |5j or 

The parameters a, r e , etc. contain information about the details of the potential 
with which the spherical wave interacts. The parameter a is known as the scattering 
length, while r e is known as the effective range. In general, an infinite number of 
terms appear in the exact expansion, but we can truncate the series for small k to 
achieve a desired accuracy with a finite number of terms. 

In a sense, we are using low-energy scattering to probe some of the properties of 
the potential. While we can use this to determine some of the details, we cannot 
determine all of them (unless we measure the phase shift exactly for all k). There 
will always be infinitely many potentials that have the same scattering length. Even 
if we did distinguish some of them by different effective ranges, there would still be 
an infinite number of potentials that share both a and r e . 

Fortunately, this is an asset and not a hindrance. What this means is that the 
low-energy behavior is insensitive to the detailed form of the potential, and any two 
potentials that give the same parameters are equally good approximations. In fact, 
sometimes it is not even necessary to find a potential to use as an approximation. 
The effective range expansion can be used to make predictions, and other derived 
quantities can be written in terms of the scattering length and other parameters. By 
experimental measurements, we can determine a and r e and then use them to predict 
other quantities. 

Let us look at an example to make some of these ideas clearer. We shall take the 
potential to be a square well 

V(r) = -V 9(l/A-r). (1.7) 
Inside the well, the wavefunction has the form 

ip(r) = A sin (Kr), (1.8) 

where 

K = yj2m(E + V Q ). (1.9) 
Outside the well, the solution takes the form 

V>(r) =Bsm(kr + 6 ), (1.10) 

where 

k = V2mE. (1.11) 
Since the quantity dhx(ip(r))/dr must be continuous, 

Kcot(K/A) = kcot(k/A + 5 ). (1.12) 

For very low energies, we can expand both sides of Eq. ()1.12j) in powers of k and then 
match terms. Notice that the phase shift occurs on one side, and Vq on the other. 
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This ensures that when we match expansions, the parameters for Sq will be written 
in terms of the parameters for the potential. 

Without going into details, we present the results of this expansion: 



1 gA 



a taa{g)-g' 

(1-aA) 2 
Te ~ a 2 A 3 ' 



1.13) 

;i.i4) 



where g = y/2mVo/ A. The one condition required to truncate the expansion is 
that J; < A, which should come as no surprise. A particle of energy k 2 /2m has 
a wavelength of about 1/k. If we demand insensitivity to the short distance behavior, 
then the wavelength must be much greater than the range of the interaction. High 
energies probe small distances, but low energies do not. The parameter A acts as a 
momentum cutoff. Incidently, this square well potential can be used to approximate 
any other short-range potential. 

By choosing g, and hence Vq, to be a function of A, we can make the scattering 
length a independent of the cutoff by always making sure that Eq. ()1.13|) is satisfied 
for the same constant a. This defines a whole set of potentials that give identical scat- 
tering lengths. The only difference between them is the value of r e . As A changes, so 
will the effective range. This is where the errors in a simple square well approximation 
come in. If our potential had another free coupling to adjust, then we could typically 
make r e independent of the cutoff, and the error moves to the next parameter. So we 
see that we can define many different potentials (one for each value of A) that can 
serve as equally accurate approximations. We can even go so far as to let A — > oo 
and end up with a delta function potential well. 

We should also take a moment to look at the relative size of our parameters. The 
scattering length appears to be around 1/A, and this would imply that r e is also 
about 1/A. [We say that a and r e are both 0(1/ A).] This is called a "natural" 
theory. All of the parameters have their lengths set by the underlying length scale 
of the problem, in this case 1/A. However, it is possible that g could be fine-tuned 
to make a much larger than 1/A. Even though aA ^> 1, the size of r e would still be 
0(1/ A). The effective range is still set by the range of the potential. However, this 
too could be made unnaturally large if there is a second parameter in our theory that 
we could fine-tune. 

In addition to looking at low-energy scattering states, we could also look at low- 
energy bound states. The condition for having a bound state of energy —B is 



g 2 A 2 - 2mB cot Jg 2 A 2 - 2mB/A = -V2mB. (1.15) 



Treating mB/A 2 as a small quantity, we expand the left-hand side to leading order 
to obtain the relation 

B = —(gAcot(g)) 2 = — f-^-r V. (1-16) 
2m vy vy;; 2m\l-aAj v ; 
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For very large a, the binding energy behaves like I /ma 2 . This shows the connection 
between large scattering lengths and shallow bound states. If we wanted to expand 
the equation for B to higher orders, we would find that for large a it is expanded in 
powers of 1/aA ~ r e /a. Once again, our expansion parameter is set by two widely 
separated length scales. 

For a natural theory, the scattering length would be set by the scale 1/A. In this 
case, Eq. implies that the binding energy is also set by the potential's length 

scale: B ~ A 2 /2m. This argument should not be taken too seriously because such a 
binding energy would violate the assumption that mB/A 2 <C 1, which was required 
to expand Eq. ()1.15|) in the first place. The leading approximation for the binding 
energy also diverges when aA = 1 showing that it breaks down for a natural theory. 
Nonetheless, the fact that B ~ A 2 /2m remains valid as can be seen by comparing 
Eqs. (HHnj) and (Tl~T5J) for the case a ~ 1/A. 

1.3 Effective Field Theory 

We now examine the interaction of two identical particles from the field theory point- 
of-view. This requires some prior knowledge of field theory Lagrangians and how 
Feynman rules are derived from them. 1 Once again, we will sacrifice rigor for sim- 
plicity and simply state many results so that readers unfamiliar with perturbative 
field theory can follow the discussion. The important things at this point are not the 
mathematical steps, but the ideas behind them. 

Suppose we have two identical particles interacting via some unknown short-range 
potential and we wish to approximate the behavior. We will start the approximation 
with the Lagrangian 

C = <jf{x) (ij t + ^ 0(f) + C [0*(x)0(f)] 2 , (1.17) 

where Cq is the coupling constant for the two-body interaction. We require that all 
interactions in the Lagrangian be local operators. That is, they are the product of 
fields at the same point. These operators can be thought of as contact interactions 
since they act only at a single point. 

If we expand the scattering amplitude perturbatively using this interaction, the 
first term will consist of a single vertex with a value of — iCq. The next term will 
contain two vertices, two propagators, and a single loop integration. Unfortunately, 
the loop integral behaves like / dk which diverges linearly. If we try calculating the 
next term, we find that it diverges quadratically. Each successive term has a more 
severe divergence than the last. It appears that our hope of a perturbative expansion 
is lost. 

^ee, for example, Ref. £Q 
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Or is it? The reason we encounter this divergence is because we have states 
with arbitrarily large differences in momenta coupled to one another. However, our 
approximation is not meant to be used for arbitrarily high momentum. It is only valid 
for low momenta that cannot resolve the internal details of the interaction. For the 
moment, let us impose a cutoff A on the momentum values. This represents the point 
at which our approximation breaks down and knowledge of the true potential becomes 
necessary. In the center-of-mass frame, the scattering amplitude as a function of the 
momentum p for the first two perturbative terms turns out to be 

A(p) = -C -^(A-mp/2). (1.18) 

To calculate A, we would need to know the values of Co and A. Yet the reason we are 
using this approximation is because we do not know the underlying potential, and we 
want our results to be insensitive to such details as the exact value of A. 
The solution is to let Co be a function of the cutoff. At p — 0, 

AO) = -Co - ^A, (1.19) 
which is related to the two-body scattering length a by 



8tt 8tt V 4tt 2 

If we invert this equation to get a perturbative expression for Co as a function of a, 
we obtain 

c^Wi-M). (1 . 21) 

This allows us to rewrite Eq. ()1.18j) as 



AP) = (1 - iap) ■ (1-22) 



m 



Our originally divergent amplitude has now been cast in terms of the finite scattering 
length, which we can determine from experiment. 

This procedure that we have just sketched is an example of renormalization. We 
regulate any divergent integrals with some sort of cutoff, allow the coupling constants 
to be functions of this cutoff, and then choose the form of the couplings to remove 
this cutoff from physical quantities. Because the coupling constants change as the 
cutoff changes, they are said to "run with the cutoff" and may be referred to as 
"running couplings." In our example, all reference to A is eliminated, at least up to 
the perturbative order in the couplings we have expanded in. In general, this will 
not be the case. Instead, there will remain some cutoff dependence involving inverse 
powers of A. The cutoff can then be taken to infinity, leaving a finite result. 

Notice that the second term in Eq. (jl.22|) contains the product ap. Our approxi- 
mation is only valid for low momentum, such that p~ l is much greater than the range 
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of the underlying interaction. If this is a natural theory, we expect that a will be of 
the same order as the interaction range and hence ap <C 1. In other words, we can 
expand our amplitude as a power series in ap, and we have seen just that. The — iap 
term is much smaller than the first term of 0(1). 

We have only included one interaction in the Lagrangian so far. With just this, we 
cannot produce a theory with a fixed, non-zero effective range. To control the effective 
range, we shall add another operator of the form C 2 V (<f)*(x)<f)(x)) ■ V (<p*(x)<p(x)). 
Since it contains derivatives of the fields, it is sometimes referred to as a "derivative 
interaction." If we were to add its contributions to A(p), we would find divergences 
that need to be regulated. This can be done once again by applying the cutoff A and 
then properly choosing (^(A). This allows us to control the effective range r e , and 
after renormalization, the divergent quantities can be rewritten in powers of r e p. For 
a natural theory, this is the same order as ap and can also be used as an expansion 
parameter. 

This example has helped to introduce the general concepts of effective field theory 
(EFT) [3ini[Zj. The main idea behind EFT is that the details of the underlying physics 
are embodied in the coupling constants of the theory. By relating the couplings to 
physical quantities (like scattering length and effective range), other calculations can 
be written in terms of these physical parameters without any knowledge of the detailed 
underlying theory. In addition, interactions containing more and more derivatives or 
fields contribute less and less to the overall result at low energies, just as higher order 
terms in the multipole expansion contribute less and less to the electrostatic potential 
far from the source. 

To see this, we must look at the dimension of each operator. We will assume that 
the range R of the underlying interaction sets the scale of the terms in our Lagrangian 
density: 

C = <J>*(x) (ij t + <f>(x) + C (0*(f)0(f)) 2 + C 2 V (0*(f)0(x)) ■ V {4>\x)4>{x)) ■ 

(1.23) 

If that is the case, then we would expect any length variables like x to scale like R, 
which we shall write as [x] = R. This implies that [d/dx] = R~ l and hence [V] = R" 1 . 
By these arguments, the kinetic energy operator would scale like [V 2 /2m] = m~ 1 R~ 2 . 
A look at our Lagrangian density shows that the time derivative must also scale in 
the same way: [d/dt] = vrC x R~' 2 . This relation then implies that the time variable 
scales as [t] = mR 2 . 

To find the scaling behavior of the fields and couplings, we observe that the action 
S = J dt J d 3 x £ is dimensionless. Since [dtd 3 x] = [tx 3 ] = mR 5 , we must have 
[£] = m _1 i?~ 5 . Therefore, every term that makes up a part of C must also scale in 
the same way: 

= m- l R~ 5 , 
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C o (0*(x)0(x)) 2 
C 2 V(0*(x)0(z)) ■ V(0*(f)0(f)) 



—5 



From this, we can easily draw the following relations 



[0] 
[Co] 
[C2] 




m- L R, 
m^R 3 . 



Since the couplings are supposed to approximate the details of the short-range 
behavior, we expect that their scale should be set by the scale of the short-range 
interaction. This is in agreement with the leading order approximation to the scat- 
tering length in Eq. (jl.20|) where a cx mCo ~ R. Notice the factor of m~ l in both Co 
and C 2 . This same factor occurs in every coupling. Every interaction operator will 
consist of fields (</>), derivatives (V), and some coupling. Every field scales as i?~ 3//2 , 
while every derivative scales like R~ x . If the fields and derivatives combined scale like 
R~ N , then the coupling must scale like m~ 1 R N ~ 5 . 

Of course, any dimensionless quantity containing the couplings must have another 
scale to balance out the powers of R. This scale comes from the momenta of the pro- 
cess being studied. If the system has a typical momentum of p, then a dimensionless 
quantity would be a function of pR. Hence, if pR C 1, we can expand in powers 
of pR. An operator of dimension R has a coupling proportional to R N ~ 5 and con- 
tributes at order (pR) N ~ 5 . This allows us to classify the operators and decide which 
Feynman diagrams are needed to achieve any given accuracy. 

Returning to our example, a single Co vertex adds a term of order R ~ a to A. 
Two Co vertices contribute R 2 p ~ a 2 p, while three add R 3 p 2 ~ a 3 p 2 . At this point, 
we notice that a single C 2 vertex is of order R 3 p 2 , implying that we must include C 2 
if we want to calculate C effects beyond the first two orders. Because the effects 
of the C 2 interaction do not contribute to leading order, we classify the operator 
V • V as irrelevant. The "irrelevant" qualifier in the term 

irrelevant operator does not carry the usual meaning, implying that this operator is 
not needed at all. As we have seen, it is vital if we desire a cutoff-independent effective 
range. Rather, we use the term irrelevant operator to mean that the operator's effects 
are smaller than the leading order effects. The renormalization group provides a more 
precise definition based on the scaling behavior of the operator about a fixed point, 
but this is more technical than required for our discussion. 

We now have a well-defined expansion method. We can calculate increasingly 
accurate results by adding additional operators, which means fitting more couplings 
to experimental data. One might wonder why this is any different from curve fitting. 
Why not create some arbitrary potential that includes a bunch of adjustable param- 
eters, and then do some sort of least-squares fit to experimental data to determine 
the optimal parameter values? The difference is that EFT provides a systematic way 
of improving the result, instead of merely guessing a potential that may or may not 
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improve things. For any given system, we can write down local operators that could 
appear as interactions, provided that they satisfy all symmetries of the system. These 
operators can then be classified based on their dimension, and the point at which they 
each contribute to the expansion is known. 

There are cases where the above expansion method is not as useful as might be 
expected. The most obvious is when a is very large. Since we require ap <C 1, a large 
scattering length severely limits the range of momentum in which our expansion might 
converge. Another case would be the study of bound states. The typical momentum 
in a bound state is of order 1/a, so once again the expansion would not converge. The 
solution in these cases is to treat ap as 0(1) and sum all terms. If this is possible, 
then we can still treat powers of r e p perturbatively for the case where r e <C a, and 
we have an increased momentum range over which the expansion is valid. 

Fortunately in our example, the contributions to A from terms containing only 
Cq form a geometric series. The result of the summation is 



S)7T 

m 



A( P ) = - 
which leads to a scattering length of 



■ill 



C n m 



2A\ 

7T / 



-, -1 



Ip 



1.24) 



/ 8vr 



2A 

7T 



\Com 

To keep the scattering length cutoff independent, we choose Co (A) to be 

2aA^ - 1 



1.25) 



C (A) 



in 



7T 



1.26) 



If we let a — > oo in the above equation, then we see that Co approaches (47r 2 )/(mA). 
It is by fine-tuning Co to be close to this value that we find a system with a very 
large scattering length. This also illustrates that the running coupling Co possesses 
a "fixed point." It is common to consider the behavior of the dimensionless versions 
of the couplings as A — ► oo. In this case, ACo — > Air 2 /m which is a fixed constant 
value. This is a typical behavior for many couplings. When we consider the three- 
body problem later, we will find that the three-body coupling exhibits another type 
of behavior: a limit cycle. Limit-cycle behavior simply means that the coupling 
approaches a periodic function as A — > oo. 



1.4 Three-Body Problem 

The EFT approach has been applied with great success to many low-energy two- 
body systems [HI EJ EE HJ- A logical next step would be the application of EFT 
techniques to the three-body problem. Recent investigations of such systems by 
Bedaque, Hammer, and van Kolck have turned up some surprising results |T2]. For 
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three identical bosons, dimensional analysis suggests that only the leading order two- 
body interaction is needed to obtain the leading order three-body results. While 
this is indeed sufficient to remove any integration divergences in the equations, cutoff 
sensitivity still remains in the form of a non- unique A — > oo limit. The remedy is 
to introduce a three-body contact interaction of higher dimension, indicating that 
short- distance effects are playing an important role. In addition, the coupling of the 
three-body interaction is found to exhibit a limit-cycle behavior, unlike the more 
common fixed-point behavior of running coupling constants. 

I choose to investigate low-energy three-body systems in a slightly different manner 
from that used by Bedaque, Hammer, and van Kolck. For sufficiently low energies, the 
potential V(r) is replaced by a series of point-like potentials, namely products of <5 3 (f) 
and its derivatives. As in the case of EFT, the coupling constants multiplying these 
terms approximate the short-range details of V(r) and must be fit to experimental 
quantities such as scattering length and effective range. 

The focus of my calculation is on the bound-state equation, which provides a com- 
plimentary approach to the scattering methods employed in other works IT2*1 [Tlj. 
My work closely follows unpublished calculations by Kenneth Wilson who developed 
the basic method and tools that I employ |T3] . I use a combination of perturbative 
expansion and numerical methods that not only provides an independent verification 
of previous results but also offers several possible advantages. There is no need to in- 
troduce a composite field, such as the "dimeron" used in Ref . > nor is it necessary 
to resum Feynman diagrams for the dimeron's propagator. I introduce a technique 
for expanding the bound state equation in inverse powers of the cutoff, giving us 
the ability to examine cutoff-dependent behavior order by order. Additionally, the 
numerical methods I use allow computation to high accuracy, typically 11 to 12 digits 
in the couplings and binding energies. 

In this dissertation, I derive the expanded set of three-body integral equations 
and present the results of my numerical calculations. I numerically "prove" that a 
three-body coupling term is necessary and sufficient to renormalize the three-body 
bound-state spectrum. The three-body coupling exhibits limit-cycle behavior, and I 
present an argument for why this should be expected. I also calculate the three-body 
bound-state spectrum and how it changes with the two-body bound-state energy. 
This data is an essential part of computing Efimov's function for three-body bound 
states to high accuracy. This function can be used to determine the binding energies 
of three-body bound states when there is a large scattering length ^H]- Efimov's 
function is currently being used to study Feshbach resonances [T7j . 

1.5 New Results 

To the best of the author's knowledge, this dissertation contains several new results, 
some of which are extensions to work done by Wilson [To] : 
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• Correction of minor errors in the leading order equations and the derivation of 
the equations for first order corrections. 

• Modification of Wilson's original leading order equations to maintain high ac- 
curacy when B 2 — B 3 in numerical calculations. 

• Detailed explanation of the derivation of the three-body integral equations and 
the use of Wilson's methods in that derivation. 

• Proof of and corrections to Wilson's assertion that certain integrals which are 
encountered when computing corrections to the leading order results are iden- 
tically zero. 

• Derivation of new integration limits to simplify numerical calculations. 

• Implementation of numerical methods to solve leading order equations. 

• Use of leading order equations to compute Efimov's universal function to high 
accuracy. 

• Use of leading order equations to analytically and numerically verify similar 
work by others. 



1.6 Organization of Dissertation 

After outlining the use of the <5 3 (r) potential in the context of the two-body problem, 
I derive some formulas in Chapter 2 that will be of use later in the analysis of the 
three-body problem. Chapter 3 gives the derivation of the three-body bound-state 
equation and outlines some assumptions made in my calculations. I digress briefly in 
Chapter 4 to explain the method used for expanding the equations perturbatively and 
then apply the technique in Chapter 5. This results in the final form of the equations 
used for numerical computation. Chapter 6 outlines some details of the numerical 
technique and discusses some of the issues that arise from the need for high accuracy. 
Following this is a discussion of analytic and numerical results in Chapter 7. Chapter 
8 is devoted to introducing Efimov's function and then numerically computing it using 
the leading order equations. I conclude in Chapter 9 with a brief discussion of areas 
for future work. 
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CHAPTER 2 



THE TWO-BODY PROBLEM 



As an introduction, we first consider the two-body problem. This provides a simple 
example that illustrates the approach we use with the three-body problem. More 
importantly though, the results derived in this chapter will carry over unaffected to 
the three-body system and are in fact needed to attack that problem. All equations 
refer to identical bosons of mass 1/2 and S-wave bound states of such particles. 
These restrictions simplify the equations while preserving the system characteristics 
we desire. 

2.1 Bound State Equation 

The general two-body Schrodinger equation is 



where Vi 2 (fi', f 2 \ f[, r^) = {f[, f 2 \ V± 2 \ f[', f^) represents the position-space matrix 
elements of the two-body interaction operator. We choose the matrix elements of V\ 2 
to be 




{fl,r 2 | V 12 \fi,r 2 ) 



g 2 8 3 {r 1 -^)5\r l > -f 2 >)5 3 



( 



-> . -> ->/,->> 
n + r 2 _n + r 2 

2 2 



) 



(2.2) 



with g 2 denoting the two-body coupling constant. Our choice ensures that 



= -g 2 j d 3 r/ S 3 (fl - r 2 )S\fl' - r 2 ')S 3 [*±* - 
= ~92 S 3 (n - f 2 ) I d 3 f[' d 3 r 2 ' 5 3 (fl' - r 2 ')6 3 ( - *-±^j ^ 




ijj(fi,f 2 ) 
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= -g2 # 3 (n - f 2 ) J d 3 f[' 5 3 y " 1 ^ r2 - rl'j n') 

= -</2 <> (n -r 2 )ip\ — - — , — - — I 

= -92 6 3 {r-l-r 2 )^{r u r 2 ), (2.3) 

which is the form we desire. In the center-of-mass frame, the bound-state equation 
for two identical bosons with binding energy — B 2 can be simplified to 

- 2 V 2 i[)(r) - g 2 <5 3 (f) i[)(r) = -B 2 i[)(f), (2.4) 
or in momentum space, 



2 P 2 <P(P) -92 J = -B^m- (2-5) 

Here we have begun the convention of labeling two-body couplings and energies with 
a "2" subscript to distinguish them from the corresponding three-body parameters 
we will encounter later, which will carry a "3" subscript. 

Equation ()2.5|) has been previously studied and is known to suffer from divergences 
for any constant, finite value of g 2 [TH]. To see this, define the constant 



Equation (|2.5|) can then be written as 

2p 2 <P(p)-g 2 l=-B 2 <P(p)- (2-7) 

This implies that 



<P{P) = d^T" (2-8) 
2p 2 + B 2 



Using Eq. O in Eq. yields 



J= / & 1 (29) 

J {2irf2q 2 + B 2 ' W 

or equivalently 

1 _ f d 3 q 1 _ dq q 2 

g 2 J (2tt) 3 2q 2 + B 2 Jo 2tt 2 2q 2 + B 2 1 ' ] 

Any finite value of B 2 results in a divergent integral, and hence a value of zero for 
g 2 . As we discussed in Chapter 1, the divergence is due to the coupling of states with 
arbitrarily large momentum differences. In EFT terms, we are including high energy 
states that are beyond the region of validity for our approximation. Therefore, it 
should come as no surprise that we encounter problems. 
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To remedy the situation, we must find a way to regulate the problem. Let us 
define 

U 2 (f) = J exp(ip • f) U 2 (p), (2.11) 

where 

U 2 (p) = exp(-p 2 /A 2 )- (2-12) 

For a finite value of p, allowing A — > oo results in U 2 (p) = 1 and hence U 2 (f) = S 3 (f). 
In essence, this defines a regulated delta function that is "smeared out" in the sense 
that it now has a range of 0(1/ A). This function can be used to regulate the two-body 
interactions by replacing the first two delta functions in Eq. (J2.2)) with U 2 , 

(n, r 2 I V 12 I rl', f 2 ') = -g 2 U 2 (f[ - f 2 ) U 2 (f[' - f 2 ') 5 3 ^l±ll _ r ±±[£^ . (2 .i 3 ) 

The momentum space matrix elements become 

= -g 2 J d 3 fl d 3 f/ exp(-ipi ■ fl) exp(i<g- ■ f/) U 2 (fi - f 2 ) 
x£/ 2 (ri -r 2 ) d f 1 



V 2 2 / 

-02 y (2tt^ 3 (270 3 U2 ^ U2 ^ rf3 ^ rf3 ^ /exp ( _ ^ 
x exp(zgj ■ rj') exp (ifc ■ (fi — f^jj exp (ik' ■ (fi' — f* 2 ] 
,3 ( n + r ~2 _ n + r2 f \ 

V 2 2 / 

d 3 £ d 3 £' 

(27r) 3 (27r) 3 ' 

x exp • fi + q 2 ■ (f[ + f 2 - f[') + k ■ (f[ - f 2 ) 

+ k ■ (2fi — fi — f 2 ) 



92 / /o nwo ^ U 2(k)U 2 (k') / dPfidPfi exp(-ipi ■ fl) 



r d 3 k d 3 P -> -> 

-92 J J^J^ U ^ U ^ k '^ 9 + <T2 + k- k') 

x5 3 (-p- 2 + q 2 -k- k') 5 3 (qi - q 2 + 2k') 

-92 J ^U 2 (k)U 2 {^^) W 5 3 (-p. + Iqi + \q 2 + k) 

x5 3 {-p 2 + |gl + \q 2 - k) 

-92 (2tt) 3 U 2 (\q 2 - \qi) U 2 (\p t - \p 2 ) 

x5 3 (qi + q 2 - fi - p 2 ), (2.14) 
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which is used to calculate 



(PilViM = f S%(^\Vu\gj)(m 



., (2tt) 6 

= -g 2 (2tt)- 3 J d 3 q 3 U 2 (\q 2 - \q[) U 2 (fpl - \p 2 ) 

x5 3 (q[ + q 2 - pi - p 2 )(j)(qj) 
= -92 U 2 (|pl - fpl) | 4%^2 (§p! + \V2 - gl) 

X0(9t,Pl +P2 - qi). (2.15) 

In the center-of-mass frame pi + p~2 = <h + <£ = so we may write the eigenvalue 
equation as 

(p\ + {-Pif) 0(Pi, -Pi) - #2 £Mp1) y U 2 (q[) (j)(qi, -q x ) = -B 2 0(p"l, -pi). 

(2.16) 

Making the substitutions pi — > p, q[ — > q, , and 0(g, — g) — > 0(g) simplifies the equation 
further, 

G? 3 (f 



2p 2 <j)(p)-g 2 U 2 (p) J ^UM<t>{q) = -B 2 (f>{p)- (2-17) 

The previous equation shows that the interaction term is of the form J dq V(p, q)<f)(q), 
where V(p, q) oc U 2 {p)U 2 (q) . Potentials that have the p and q dependence in separate 
factors are called "separable potentials." The separable nature of this potential will 
be of use later, and this is why we have chosen this particular form for the regulated 
two-body interaction. 1 It is important to note that how we regulate high-momentum 
states should not matter. This allows us to choose a regulator that is convenient 
for calculations. Just as the details of the short-range interaction are unknown and 
unimportant, details of how we regulate the short-range interaction are not important 
for low-energy physics. 

To find the new expression for l/g 2 , we need only alter our previous arguments 
slightly. Let 

J(A) = / ^y UMm ' (2 ' 18) 

so that Eq. ()2.17j) becomes 

2p 2 0(p) -g 2 U 2 (p)l(A) = -B 2 <j)(p), (2.19) 

leading to 



1 For an early example of the use of separable potentials in nuclear physics, see Refs. ^H] and [23j- 
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Once again, using the above form of in the definition of 1(A) leads to the formula 
for g 2 : 



U *® % (2.21) 



g 2 (A) J {2Tif2q 2 + B 2 

d 3 q exp(-2g 2 /A 2 



(2vr) 3 2q 2 + B 2 



(2.22) 



The gaussian factor now renders the integration finite. By choosing g 2 to satisfy this 
equation for any value of the cutoff, the bound-state energy B 2 is guaranteed to be 
independent of A. For very large values of A, we may expand g 2 as a power series in 

8V2tt 3 / 2 f_ ^jB~ 2 x 



92 = -^j— ^1 + + ..■)■ (2.23) 

The dimensionless two-body coupling G 2 = A g 2 exhibits fixed-point behavior, ap- 
proaching 8v^2tt 3//2 as A — > oo. The fact that this constant is not small is an indi- 
cation that the problem is nonperturbative and an expansion in powers of g 2 is not 
feasible. 



2.2 Irrelevant Operators 

As an EFT, the current theory is capable of removing all cutoff dependence from one 
physical observable which we have chosen to be the two-body bound-state energy. 
However, any other observable quantity will contain cutoff dependence, and the accu- 
racy can be no better that 0{\/B~ 2 / A). 2 Improving the accuracy and removing cutoff 
dependence from a second observable requires the introduction of other operators in 
the effective potential. 

To demonstrate this, we attempt to remove the cutoff dependence from two pa- 
rameters: the two-body scattering length a and effective range r e . We can introduce 
a second operator by redefining U 2 (p) to be 

U 2 (p) = (l + h 2 ^j exp(-p 2 /A 2 ). (2.24) 

We have chosen this form for convenience, but EFT states that the short-distance 
details are unimportant. Another form could work equally well. The steps leading 
to Eq. (|2.21|) in the previous section remain unchanged, the only difference being a 
dependence on h 2 in the final result. 

Roughly speaking, the introduction of this new term adds an irrelevant operator 
of the form V 2 5 3 (r) to the effective potential with a coupling strength of g 2 h 2 . This 

2 We will exhibit this shortly with the two-body effective range. 
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operator is labelled as irrelevant because we will see that the value of h 2 has no effect 
on the leading order behavior of the two- or three-body system. 

Scattering properties for two identical particles of mass m (21 El can be derived 
from the T matrix, which is a solution to the Lippmann-Schwinger equation: 

rpf , ^ \r( i \ j_ [ d ^ V(p',q)T(q,p;E) 

T(p,p;E) = V{p,p)+ — — — — . 2.25 

The +ie prescription implies boundary conditions for which there are no incoming 
spherical waves and is needed to avoid poles in the integrand. If we use a principal 
value prescription instead, which corresponds to incoming and outgoing spherical 
waves, the solution would be the K matrix: 

, , f d 3 q V{p',q)K{q,p;E) 

K(p,p;E) = V(p,p)+Vj — - _ q2/m . (2.26) 

The K matrix contains the same scattering information as the T matrix. The reason 
for using the K matrix is its simple relation to the effective range expansion 

m , m ( 1 1 9 \ , 

--pcot(5o( P )) = -- (-- + -r eP 2 + ...), (2.27) 



K(p, p;p 2 /m) 

where 5q(p) is the S-wave scattering phase shift. Because we have a separable potential 
in momentum space, we can easily obtain an explicit expression for the K matrix: 

tf(p>;£) = -^g!p, (2 . 28) 

where 



To relate the couplings g 2 and h 2 to the two-body parameters a and r e , we must 
expand K(p,p;p 2 /m) and hence D{p 2 /m) to the proper order in p. Substituting 
m = 1/2 and writing the expansion of D(2p 2 ) as V + T>\p 2 leads to the following 
expansion for 1/K: 



K(p,p;2p 2 ) u V A 2 

with 



= 1 _A(3^ + 8A, + 16) 3 

92 128^/2773/2 ' V ^ 

h\ + 8h 2 - 16 . . 

V\ = — - — jJ-, . 2.32 

32V2vr 3 / 2 A V ; 
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After equating Eqs. (|2.27|) and f)2.30|) . some algebraic manipulation yields 



1 1 

g 2 8im 



A(3/i| + 8h 2 + 16) 



128V2tt 3 / 2 



0. 



1 



2ttA ah 2 



h 2 + 



ah 2 



2ttA 



0. 



(2.33) 



(2.34) 



These two equations allow us to choose the cutoff dependent forms of g 2 and h 2 
in such a way that a and r e are cutoff independent and match the experimental 
data. In addition, the error in other A-dependent observable quantities will now be of 
0(B 2 /A 2 ). These equations also allow us to see what would happen if h 2 were equal 
to zero. In that case, the effective range would be 



ah 2 



2ttA 



1 - 



2aA 



(2.35) 



showing that it now possesses cutoff dependence and vanishes as A — ► oo. It also 
contains sub-leading corrections of 0{l/aK) ~ 0(y/B 2 '/A). 

In order to examine the limiting behavior of the couplings, we must keep in mind 
one vital relation between r e and A. For any finite value of r e , at very large cutoffs 
Eq. (jHHg) looks like 

h 2 2 + 8h 2 + 2v / 2TAr e - 16 = 0. (2.36) 



There are now three cases to consider: r e 
If r e = 0, then Eq. ()2.36|) reduces to 



0, r e > 0, and r e < 0. 



hi + 8h 2 - 16 



0. 



(2.37) 



The solution is h 2 = 4(— 1 ± Vz), indicating that h 2 is a fixed point. The equation 
for g 2 has changed, but its qualitative behavior has not. There remains a fixed point 
for the dimensionless coupling, but at a slightly different value: 



G 2 



16 + 8h* 2 + 3h* 2 



(2.38) 



Here, h\ is the fixed point value of h 2 . 

If r e > 0, the solution to Eq. ()2.36|) for h 2 will only be real if A < 16/(V^7rr e ). 
This implies that the cutoff has a limit on its value, and we cannot take the limit 
A — > oo. We must always have a finite value for the cutoff. 

Finally, for the r e < case, there will always be a solution for h 2 : 



h 2 = -A± V 32 - 2V27rr P A 



(2.39) 



As A — » oo, the solution for h 2 diverges like y/X. From Eq. 1)2.33)1 . we see that this 
behavior implies that g 2 ^ 0. 
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If we had the desire to fix any more parameters in the effect range expansion, 
we would need more tunable couplings in our interaction. This is easily done by 
extending the polynomial in Eq. (J2.24j) : 

Ui(p) = (l + h 2 ^ 2 + h> 2 ^ + ■■) exp(- P 2 /A 2 ). (2.40) 

Each additional term also reduces the errors in the other cutoff-dependent quantities, 
and there is a one-to-one correspondence between the terms in the expansion of U2 
and the terms in the effective range expansion in Eq. 1)2.27)1 . 
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CHAPTER 3 



THE THREE-BODY EQUATION 



The previous chapter laid the foundation for a two-particle system. We are now 
in a position to investigate the behavior of the system when a third identical boson 
is added. The interaction between these particles will consist of all pair-wise two- 
body interactions, as well as a three-body contact interaction that takes the form of 
a product of two delta functions. 

We begin with the position-space representation of the Hamiltonian. The interac- 
tions are regulated by a cutoff A to remove divergent behavior, with the three-body 
regulated interaction being modeled after the two-body version. The bound-state 
equation is then converted to momentum space, where it is rewritten in the form of 
a one-dimensional integral equation. At the end, we briefly discuss the necessity of 
including a three-body contact interaction. 

3.1 Position Space Representation 

In position space, the three-body bound-state equation takes the form 

- B^{f[,i%,i%) = [-V?-V|-Vl 

- 92 5 3 (n - r 2 ) - g 2 5 3 (f* 2 - r 3 ) - g 2 5 3 (r^ - f[) 

+ 93^ 3 (n -r 2 )S 3 (r 2 -r z )\ $(f[,r2,i%). (3.1) 

Here, B 3 is the three-body bound-state energy, and g% is the coupling strength of 
a three-body contact interaction which acts only when all three particles are at the 
same point. The form of the three-body interaction is somewhat arbitrary. We have 
chosen the product of two delta functions because it is the simplest one that is non- 
zero only when f\ = f 2 = 7%. Other forms could be used, but the results will remain 
unchanged due to the principles of EFT. 

In order to facilitate the transition to momentum space, we present a more de- 
tailed look at the position-space matrix elements of all interactions. All two-body 
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interactions are exactly the same as the ones used in Chapter |21 the only difference 
being the presence of a third position vector. For example, the operator Vi 2 now has 
matrix elements 

(rl, rl, r% \ V 12 \ n, r 2 , f 3 ) = (rl, rl | V 12 \ n, r 2 ) 5 3 (f 3 - r 3 '), (3.2) 

where (rl, f 2 \ V\ 2 \ rl', f 2 ) is defined in Eq. (12. 2j) . The addition of a third particle adds 
an extra delta function. The expressions for V23 and V 3 \ are obtained by a simple 
permutation of indices. Because the two-body interaction remains unchanged, so does 
the renormalization of g 2 . Its behavior is fixed by the two-body sector and carried 
over directly to the three-body sector. 

The three-body contact interaction, V123, will have matrix elements of the form 

(rl I ^123 I r/) = g 3 J d 3 f c f[ 5 3 (f t - f c )5 3 (f/ - f c ). (3.3) 

This is chosen so that 

(r t I ^123 |V> = / d^'in I V 123 I r-/)(r/\tfj) 



g 3 J d 3 f- d 3 f c 



3 

Yl 5 3 {fi-f c )5 3 {f/-f c 



^(n,r 2 ,i%') 



= 93 j d 3 f c 8 3 {f[ - r c )5 3 (f 2 - f c )5 3 (f 3 - f c )ip(f c , f ci f c ) 

= 93 $ 3 (ri ~ r 2 )5 3 (f 2 - f^)ip{f[, r 2 , r 3 ). (3.4) 

We know from the previous chapter that this potential must be regulated to isolate 
and remove the spurious short-range features that cause divergences. This is done in 
a manner similar to the two-body case. We define a new set of functions 

U 3 (p) = (l + /^exp(V/A 2 ), (3.5) 

U 3 {r) = J -!0^ex.p(ip-r)U 3 (p), (3.6) 

to regulate the three-body interaction. The regulated matrix elements of V123 become 

r 3 _ 

(n I Vi23 I f/) = g 3 J d 3 f c 1] U 3 (n - f c )U 3 (r/ - f c ). (3.7) 

The parameter h 3 plays a role analogous to h 2 , allowing the addition of a second three- 
body operator. We can use h 3 to tune another three-body observable in much the 
same way that h 2 tunes the two-body effective range. In both cases, the additional 
contribution is suppressed like an irrelevant operator. Because of this, we will see 
that h 3 can be set to zero in the leading order calculation. It is included here for 
completeness. 
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3.2 Momentum Space Representation 

Knowing the complete regulated position-space representation of the Hamiltonian, 
we can derive the momentum-space Schrodinger equation that serves as the basis for 
all of our calculations. The left-hand side of Eq. ()3.1|) and the kinetic terms on the 
right side are easily converted: 

- B 3 ip(f u f 2 ,f 3 ) -> -S 3 0(pl,pl,pl), (3.8) 
(-V?-V!-V!)^(r!,7S,rS) -> (pi+p 2 2 +pl)(f)(fi,P2,P3). (3.9) 

For the two-body interaction V±2, the momentum-space matrix elements are 



= -02 J d 3 fi d 3 f*/ exp(-ipi ■ fi) exp(iqj ■ f-)tj 2 {f 1 - r 2 ) 
x£/ a (ri -r 2 )5 I — Id (r 3 - r 3 ) 

x exp(igj • rj') exp uA; • (rj — r^n exp • (r"|' — r-j')) 
( — g g ) 5 ( r3_r3 ) 

= y ^f{2nf U2 ^ U2 ^ I ^ d3r "i' ex p(-^ • ^ 

x exp [q{ ■ fx + q 2 ■ (f[ + r 2 - f{) + q 3 ■ f 3 + k ■ (f[ - f 2 ) 
+ k' ■ (2rT — j^ — f 2 ) 



-92 J (L)3 ( 2 7r) 3 ^(fc)^(fc , )(27T) 12 ^(-pl + g + k ~ k') 

x5 3 (-p 2 + q 2 -k- k')5 3 (q-[ - q 2 + 2k')5 3 (-p 3 + g) 

-fti / (Ss^tfa - fi/2) (27r) 9 5 3 (-p 1 + §tf + §<g + fc) 

x5 3 (-p 2 + §gl + \(h - k)5 3 (-p~3 + g) 

-92 (2tt) 6 E7 2 - \qi) U 2 (\p[ - \p 2 ) 5 3 (-p 3 + q 3 ) 

x5 3 (qi + q 2 - fx - f 2 ), (3.10) 



which is used to calculate 



(2tt 
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-g 2 (27r)- 3 | d 3 q 3 U 2 (\q 2 - \q[) U 2 - \p 2 
xS 3 (-p 3 + q 3 )5 3 (q{ +q 2 - p[ - p 2 )<p(qj) 



= -92 U 2 - \p 2 ) J -^r 3 U 2 (§pl + \p 2 - q x 

X0(9l,Pl+P2-9l,P3)- (3-11) 

Because the bound state consists of identical bosons, in the center-of-mass frame we 
know that p[ + p 2 + p 3 = and <f>(qi,p~i + p 2 - qi,Pz) = 4>(ql,p~3,P~i +Pi ~ <fi)- This 
allows us to write 

(Pi\Vi 2 \4>) = -g 2 U 2 (|pl - \p 2 } J jr^ u 2 (ql + |ps) 0(gi,P3, -qi ~ Pz)- (3.12) 

For reasons that will become apparent shortly, we define the following "pseudo- 
wavefunction" : 

$(p) = / 7^33 t W+ lp)<t>{£p,-q-p)- (3.13) 

Equation (|3.12|) is now written as 

(PilK 2 |0) = -<& 1/ 2 - ipi) $ (^), (3.i4) 



with the other two-body interaction terms being obtained by a cyclic permutation of 
indices: 

(pllW) = -g 2 U 2 (±p 2 - ipl) $(pl) = -<? 2 C/ 2 (pi + ipi) $(pl) (3.15) 

(p1|V 3 i|0> = -£ 2 f/a (|P3 - \fi) HP2) = -92 U 2 (fx + \p 2 ) $(p 2 ) (3.16) 
Turning to the three-body interaction, the matrix elements are 

<*i«.«> - j^r;<m><f t ^>< m > 

= 93 J d 3 fl d 3 f- exp(-ipi ■ fl) exp(ig* • r/) 
x / d 3 f c n U 3 (fi - f c ) Uzifj - f c ) 
f d 3 k ■ d 3 k • f 

= 93 J (2tt)9 (2tt) 9 J d3fi ^ d3r " cex PH^ ' ^ e Mi<lj ■ r/) 

x exp (iki ■ (fl - f c )j exp fikj ■ (f/ - r c )\ U 3 (ki)U 3 (kj ) 
= g 3 [ d% d^k-U^U^k-^f 5 3 (-pi + k x )5\-p 2 + k 2 ) 



x5 3 (-p 3 + k 3 )5 3 (qi + h )6 3 (q 2 + k 2 )6 3 (q 3 + k 3 ) 
x5 3 (kx + k 2 + k 3 + kx + k 2 + k 3 ) 

g 3 (2tt) 3 C/ 3 (Pi)t/3(gl)5 3 (pi + pi + p 3 - <fi ~ €2 - (3.17) 
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This implies that 



xU 3 (l2Pi ~ ft ~ ^1 ^ fft' ft; XjP* - ft - ft) 

= gMp{)u 3 (p 2 )u 3 (f 3 ) [ fy^fJlM &M) 



(2tt) 3 (2tt) 3 

xU 3 (-q[ - gl)0(gl, ^, -gl - gl), (3.18) 



where the last step uses the fact that Y$=iPi = 0. If we define 

(2-nf (2tt) 3 (ft ) ^ 3 (~ft _ ft)<Kft> ft' ~ft ~ ft)' ( 3 - 19 ) 
then Eq. ()3.18|) may be written simply as 

(muM = gsUMWz&WM)®!- (3.20) 

Overall, the momentum-space three-body bound-state equation takes the form 

- B 3 (f)(pup 2} p 3 ) = (p? +pj+ pl)0(pi,P2,P3) - 92 U 2 (f 2 + $(pi) 

- #2 E/a (pi + §pl) ^(Pl) - ft t/ 2 (fpi - |pl) $(p1) 

+ gsUsipDU^UM)^. (3.21) 

3.3 Elimination of Wavefunction 

The bound-state equation can be simplified further by choosing to work entirely with 
the pseudo-wavefunction, $(p), rather than the wavefunction 4>(pi,p~2,P~3)- Because 
we limit our investigation to S-wave bound states, $(p) is only a function of the mo- 
mentum's magnitude, not direction. The result is a one-dimensional integral equation, 
as opposed to the multi-dimensional equation required for 4>(p[,p~2,P~3)- 

To eliminate the wavefunction, we solve ()3.21|) for (f> and substitute the result into 
(EH: 

$ p = / 7^3 2^ 2,^,^2, p ^2^2(p+|g)^(g 

7 (27r) 3 p 2 + g 2 + (p + q) 2 + B 3 L V ^/ 

+ ft U 2 (q + \p) $(p) + ft 27 2 (|g - ip) $(-p - g) 

- ft U 3 {q)U 3 {p)U 3 {-p- g)$x] , (3.22) 
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which can be rewritten as 

1-02 



'+ W) 2 



( 2n )z p 2 + q 2 + (p + q)2 + B 3 \ 

d 3 q U 2 {q+\p)U 2 (p 



92 



(2tt) 3 p 2 + g 2 + (p + g) 2 + B, 



$(g) 



+ 92 



~93 



d 3 q U 2 {q + \p)U2{\q~ \p) , - ^ 
(2tt) 3 p 2 + g 2 + (p + g ) 2 + 5 3 1 P ? ' 

rfV u 2 ti+y?)u 3 ®u 3 (p)u 3 (-f-?] 



(27T) 



P 2 + g 2 + (p + g) 2 + S 3 



$i. (3.23) 



By performing the change of variables q^q—^P, the left-hand side of the previous 
equation becomes 



I -92 J 



(2vr) 3 2g 2 



\v 2 



B 3 



®(p)=g 2 D(-lp 2 -B 3 )<S>(p), 



(3.24) 



where D(E) is defined in Eq. (I2~2HJ) . 

Next, consider the second integral on the right-hand side of Eq. ()3.23|) . The change 
of variables q — > — g '— p allows us to rewrite the integral as 



C/2 



d 3 g- C/ 2 (-g- |p)^ 2 (-|g 



P. 



(2tt) 3 2p 2 + 2g 2 + 2p • g + £3 



7 (27r)' ;i p 2 + g 2 + (p + g) 2 + 5 3 

(3.25) 

which is exactly the same form as the first integral on the right-hand side of the same 
equation. The two can therefore be combined to simplify the equation. 
Finally, we define 



Di(p) = 



d 3 q U 2 (q+ W)Us(q)U3(p)U 3 (-p 



Q 



(27T) 



2p 2 + 2g 2 + 2p-q + B 3 



(3.26) 



so that the three-body interaction in Eq. ()3.23|) can be written as — g 3 D 1 (p)Q 1 . 
Dropping any unnecessary vector notation, the bound-state equation takes its final 
form 



$(p) 



D 



-\p 2 -B 3 
g 3 D 1 {p)<& l 



q 2 dq 

~^ 2 J -i ~~ 2p 2 + 2q 2 + 2pqz + B 3 



92 D 



-IP 2 - B; 



(3.27) 



where p-q — pqz. The definitions for quantities in the above equation are reproduced 
below: 



$(p) 



/d?q 
j—^U 2 (q+^p)(f)(q,p,-q-p), 



(2*) 



(3.28) 
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D{E) 



1 

92 

1 



-V 



d 3 q U 2 {qf 



poo rl 

4^^ jo J-i 



(2tt) 3 2q 2 -E' 

U2(q + \p)UMUz(p)Uz(p + q\ 



2p 2 + 2q 2 + 2pqz + B 3 



d 3 q d 3 p 
(2tt) 3 (2tt) 3 



U 3 (q)U 3 (p)U 3 (q+p)<f)(q,p,-q-p] 



(3.29) 
(3.30) 
(3.31) 



At this point, the only thing preventing the exclusive use of $(p) in our calculations 
is the fact that the definition of $i still involves <fi. To remove this dependence, we 
insert the value of obtained from Eq. (|3.21J) into (|3.31jl : 



d 3 p U 3 {q)U 3 {p)U 3 {q + p\ 



(2tt) 3 (2tt) 3 2p 2 + 2q 2 + 2pqz + B 3 



+92 U 2 (q + \p)$(p) + g 2 U 2 {\q- \p)${-p - £ 
-03 U3{q)U 3 (p)U 3 {q + p )$! 



(3.32) 



Upon defining 
£ 2 



*tt 4 Jo PdP Jo Q Q l-i d 2p 2 + 2q 2 + 2pqz + S 3 ' ( ' 



the ^3 term reduces to —g 3 D 2 $i. The three #2 terms in ()3.32|) can all be written in 
the same form, provided that we swap p and q in the second integral, and make the 
variable substitution p — ► — p — <f in the third. This common form is 



92 



d 3 g $ , . /• rf 3 p ^ 3 (g)^3(p)^3(g + p)^ 2 (p+fg ; 



(2vr 



(2vr) : 



2^2 _|_ 2 g2 _|_ 2 pg Z _|_ ^ 



02 

2vr 2 



The equation for $i can now be cast entirely in terms of 



3^2 

2tt 2 7o 



rfg q 2 D l {q)<&{q) - g 3 D 2 <& 1 . 



dqq 2 Di(g)$(g). 

(3.34) 



(3.35) 



Note that the quantities $(p), $i, -Di(p), and _D 2 contain dependence on £> 2 , £>3, 
and/or A even though it is not explicitly indicated. 

For the reader's convenience, the three-body bound state equation is given below 
along with any necessary supplementary equations. 
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$(p) 



D(E) 

Di(p) 
$1 

^2 



D 



Ip 2 



92 



1 

92 



-V 



ip 2 



q 2 dq 
4^ 



U2(q) 2 



1 dz U2(q+W) U 2(P + |g) 
-l ^ 2p 2 + 2g 2 + 2pgz + 5 3 



(2vr) 3 2g 2 - £ 



471-2 
3^2 

2tt 2 
1 

8^ Jo 



q 2 dq 



1 d + \p)uMUz(p)Uz{p + g) 

-i * 2p 2 + 2g 2 + 2pqz + 5 3 



p 2 dp / g 2 <ig 



1 dz UM 2 uM 2 Uz{q + p) 2 

-i ^ 2p 2 + 2g 2 + 2pqz + fi 3 



3.4 Necessity of Three-Body Interaction 

Up to this point, we have implicitly assumed that a three-body contact interaction 
is a necessary part of the bound-state equation without justifying why this may 
or may not be true. In the absence of g 3 , the bound-state equation involves only 
the quantities -B 3 , and A. 1 If the two-body equation has been renormalized by 
fixing B 2 , then the only adjustable parameters are _B 3 and A. For any fixed value 
of A, the three-body bound-state equation becomes an eigenvalue equation for B 3 . 
This equation determines the entire discrete bound-state spectrum. Fitting B 3 to 
experimental data without g 3 would require fine-tuning A to an appropriate value. 
The three-body bound-state energy would be inherently tied to the value of A. For 
the additional case of B 2 = 0, the equation requires that B 3 /A 2 be a constant, causing 
the three-body bound-state energies to diverge as A 2 . 2 

To remove cutoff dependence, we introduce g 3 . By allowing g 3 to be a function 
of A, any change in the cutoff can be compensated by a change in the three-body 
coupling, keeping the bound-state energy B 3 fixed. Notice, however, that only one 
bound-state energy can be made cutoff independent in this manner. Other binding 
energies in the spectrum will contain A dependence. Our ultimate goal is to study this 
dependence, and in the next chapter we develop a perturbative expansion technique 
for doing so. 



1 It also contains <?2, which is constrained by the two-body equation to be a function of Bi and A. 

2 Taking the limit A — > oo results in a three-body spectrum that is unbounded from below. This 
is exactly the behavior Thomas found for zero-range forces |21| . 
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CHAPTER 4 



PERTURB ATIVE UNIFORM EXPANSION 

METHOD 



Our goal is to expand the three-body bound-state equation about an infinite cutoff 
and then study finite cutoff effects perturbatively. This first requires a zeroth order 
solution to Eq. (|H.27jl . We assume that B2 and B3 are much less than A 2 , making an 
expansion in powers of y/B^/A and y/B^/A possible. However, the momentum covers 
the entire range of values from to infinity If p ~ \/B~l <C A, we can obviously use 
p/A as another expansion parameter. But how should we handle the expansion when 
\fB3~ < p < A or \fB~j, <C p ~ A? And how can we expand a function that will end 
up being integrated over all three regions? 

The solution is to separately expand the function in each of the above re- 
gions. Together, the three approximations can be combined into one function that 
is uniformly valid. This means that for any value of p the approximation will be 
accurate to some chosen order in \fB~2l A or \/B~l/ A. Although it may appear we 
are complicating matters by trading one integral equation for three coupled integral 
equations, we will see that to leading order the set of integral equations becomes 
uncoupled. One of the equations can be solved analytically, while the other two must 
be solved numerically. Because the expanded equations are simpler than the original, 
we can obtain highly accurate numerical results. The high accuracy is necessary to 
ensure that any corrections to the leading order equations do not become obscured 
by numerical errors. For certain cases of atomic scattering, this could be O(10~ 6 ), 
or even smaller if there is an external magnetic field tuned to a Feshbach resonance 

Towards this end, we develop a method for finding the perturbative uniform ex- 
pansion of any function containing two widely separated energy scales, r] and A, such 
that ?; < A. 1 The following sections provide a general prescription for carrying out 

lr Thc procedure outlined here was provided by Kenneth Wilson. It is unknown to the author 
whether this represents original work. 
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this expansion to any order in rj/A. It is applicable to functions involving one or two 
additional momentum variables, as well as integrals of such functions. 



4.1 Three Parameter Function 

Consider a function of three variables, f(r],p,A), where rj <C A. We wish to expand 
the function so that, for any given value of p, the approximation is equal to / to 
within some given order of rj/A. In addition, the integral / dp f(r),p,A) should be 
accurate to the same order. 

To achieve this, we expand f(r],p,A) in three different regions. The first region 
is p ~ rj <C A, where we expand in powers of rj/A and p/A. We call this expansion 
fi(rj,p,A) using the "/" subscript to designate the "low-momentum" region. In the 
region i] <C A ~ p, the expansion parameters are rj/A and rj/p. We refer to this 
expansion as fh{f],p,A) using h to indicate the "high-momentum" range. The final 
region is rj <C p <C A. Here we expand in powers of r]/p, r]/A, and p/A. The resulting 
function will be called fd(r],p,A), where d stands for the "double expansion" region. 
The approximation to the original function can then be written as 



fi(v,P, A) + f h (r],p,A) - fd(v,P, A) ~ f(r],p, A). 



(4.1) 



We do not offer a formal proof of this statement. Instead, we illustrate the concepts 
with an example, and present arguments to convince the reader that the method is 
generally applicable. 

Suppose we desire an approximation to the function 



f(v,P,A) = 



(p + rj)(p + A) 



(4.2) 



that is accurate to 0(rf/ A 2 ). To calculate all terms of order up to and including 
i] 2 /A 2 , rjp/A, and p 2 / A 2 are kept: 



fi(v,P,A) = 



A(p + rj) 



1 _ E. + EL 
A A 2 



(4.3) 



For f h , we must keep all terms smaller than or equal to orders r] 2 /A 2 , r] 2 /p 2 , and 
rf/pA: 

f h (r],p,A) 



, _ V rf 

P p A 



(4.4) 

pijp + A) [ p p 2 

To find the double expansion, fa, we expand fi by treating p as though it were of the 
same magnitude as A. The result is 



f d (r],p,A) = - 



1 _P,P_ 
A A 2 



1 _ V T_ 

9 

P p Z 



(4.5) 
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Equally valid is expanding f h by treating p as if it were of the same magnitude as T], 
leading to 



fd(v,P, A ) 



1 

P 



r L + r L 



1 

A 



1 _ P_ ,P_ 
A A 2 



(4.6) 



p p z 

Notice that both approaches yield the same result, which should always be true. 

To prove that the approximation is accurate to the order stated, we calculate the 
difference between the function and its approximation: 



f — fi — fh + fd 



7}- 



A 3 (p + r])(p + A) A 3 
The error is smaller than the function / by a factor of rf / A 3 for any value of p. 



(4.7) 



4.2 Integrating a Three Parameter Function 

The integral J °° dp f(t],p,A) will result in a function of r] and A. Integrating 
Jo°° dp (fi + fh — fd) should result in the same function to within the accuracy of 
the approximation itself, which in our example is 0(r] 2 /A 2 ). The benefit of integrat- 
ing each term individually is simplicity. Obtaining a closed form solution for the 
integral of / may prove difficult, if not impossible. If one is obtained, then it must 
be expanded to the appropriate order. On the other hand, fi, fh, and fd may have 
much simpler forms than the original function, making them easier to integrate. After 
integration, no additional expansion is necessary since each term is already expanded 
to the proper order. 

Unfortunately, a problem arises when integrating individual terms. When inte- 
grating fi from Eq. (|4.3|) . we encounter the term 

dp A^y (4 - 8) 

which diverges as A — > oo. This quadratic divergence is canceled by an equal and 
opposite quadratic divergence from the p/A 3 term in fd- Even so, the resulting 
difference of these terms is — (rjp)/(A 3 (p + i])) which now produces a linear divergence. 
This, in turn, is canceled by the —r]/A 3 term from fd, resulting in a difference of 
r] 2 /(A 3 (p + if)) that is logarithmically divergent. Once again, a term from fd, namely 
rj 2 /(pA 3 ), cancels this divergence leaving — rf/{pA 3 (p + if)). This final term is finite 
when integrated to infinity. Unfortunately, it is now afflicted with a divergence as 
p —>■ 0. This divergence is again canceled (along with all other ultraviolet and infrared 
divergences) since the integral of the sum is assumed to converge. 

In anticipation of the fact that all divergences must cancel, we define a standard- 
ized subtraction scheme that can be applied uniformly to every integration. This will 
render all integrals finite by effectively adding and subtracting the same divergent 
term to different parts of the expansion. The subtraction scheme we utilize is defined 
as follows: 
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Integrals of positive powers of p and constant terms independent of p are 
dropped. 

A 1/p integration that diverges as p — > oo will have an equal 1/p integration 
from A to oo subtracted out. 



These rules handle all p — > oo divergences. For divergences resulting from p 
previous pair of rules is augmented with another similar pair: 

• Integrals of negative powers of p (other than 1/p) are dropped. 



0, the 



A 1/p integration that diverges as p 
from to r] subtracted out. 



will have an equal 1/p integration 



If the integral / dp f(rj,p, A) is finite, then any divergence occurring in one of the 
terms of / dp (fi + fh — fd) must be canceled by an equal and opposite divergence 
in another term to guarantee a finite result. By applying the subtraction scheme 
to the first divergence, we render it finite. However, the same subtraction scheme 
would also need to be applied to its "partner" divergence resulting in an equal and 
opposite subtraction term. The net result is that the sum of all terms before and 
after the subtractions are applied must be equal. This same argument could be 
used to formulate alternate schemes that would work just as well, provided that they 
are applied uniformly to all integrations. Of course, formulating and applying any 
subtraction scheme requires some analytic knowledge of the integrand. 
As an illustration of these rules, an integral like 



o 



dp 



P 



A 3 (p + T])' 



which diverges as p 



oo would now become 

n2 



dp 



P _ P_ + JL 
A 3 (p + i]) A 3 A 3 



Ap^O divergent integral such as 



dp 



o p 2 (p + A) 



would be modified to 



dp 



V 



p 2 (p + A) p 2 A 



dp 



pA 3 ' 



+ 



" dp — 
o pA 2 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



We can also apply the expansion method to the case of gi- From Eq. (j2.21jl we 
have 

1 r d 3 q U 2 (q) 2 



92 



(2vr) 3 2q 2 + B 2 



(4.13) 
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Not only will this example help clarify the ideas presented, but it is also easily com- 
pared to the previously derived expansion for g 2 . We define 



p 2 U 2 {pf 



2tt 2 (2p 2 + r] 2 ) 



(4.14) 



so that l/<72 = Jo° dp f(rj,p, A) with 77 = \[B~ 2 . We desire an answer accurate to 
0(r]/A). The values for f h f h , and fa are: 



fi(v,P, A ) 
f h (r],p,A) 
fd(v,P,A) 



2tt 2 (2p 2 + r] 2 ) 



471-2 

1 

47T 2 ' 



U 2 (pf 



(4.15) 

(4.16) 
(4.17) 



Next, we integrate each of the above functions, applying the subtraction scheme as 
needed: 



dp 



roo 

/ dp fi(ri,p,A) - 
Jo 

roo roo | 

J Q dp f h (77, p, A) -> y dp — If 2 (p) 



2^2 _|_ ^2 



8^ 



7T 



dp fd(v,P, A ) 







dp 



A(3/if, + 8/i 2 + 16) 

128^2 7T 3 /2 

1 1 



4tt 2 



4tt 5 



0. 



The final result is 



or upon inverting, 



92 



1 A{3h 2 + 8h 2 + 16) 
lh ~~ 128^ 7T 3 / 2 

128V/2 7T 3 / 2 



A(3h4 + 8h 2 + 16) 



16-y/7T77 



A(3/i| + 8/i 2 + 16) 



(4.18) 

(4.19) 
(4.20) 

(4.21) 
(4.22) 



This matches exactly with Eq. ()2.23|) when h 2 = 0. Notice that the subtraction 
scheme terms in Eqs. f!4.18j) and (J4.20)) are equal and would cancel out if all terms 
were combined. The simplicity of this approach is also evident considering that the 
alternative would mean integrating f(i],p,A) to obtain (for h 2 = 0) 



8V2 7T 3 / 2 



A — 7]y/lTexp 



A 2 



f]\^T exp 



A 2 



Erf 



(4.23) 



and then expanding this quantity in powers of 77/ A. A non-zero value of h 2 would be 
more complicated still. 

There is one important point to note. When applying this expansion technique, it 
is possible that two terms of a given order can be multiplied to yield a term of the same 
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Label 


Sign 


Region 


Parameters 


fl 


+ 


(V,P) <A 


r)/A,p/A 


fd 




77 < p < A 


rj/p,r}/A,p/A 


fh 


+ 


77 « (p,A) 


T]/p,7]/A 



Table 4.1: Expansion terms for a function of three variables. 



order. For example, when expanding to order 77/A in the region t] ^ p ^ A, products 
of T)/p and p/A may be obtained. Each term individually must be kept because it is 
considered 0(r]/A). However, the product (rj/p)(p/A) is also 0(rj/A) and cannot be 
dropped. This is an important difference compared to most perturbative expansions 
where the product of two first order terms must be of second order. Fortunately, this 
caveat only needs to be considered in the mid-momentum region, 77 <C p <C A. 

4.3 Four Parameter Function 

To expand our three-body equation, we need a more complex version of the above 
procedure that can be applied to functions of four variables. We desire a uniform 
expansion of the function u(rj,p, q, A), where 77/A <C 1 and p and q are arbitrary. The 
steps involved are very similar to the case of three variables: expand the function 
in several regions and combine the results. Unlike the three variable case where 
we needed to expand in only three regions, we must now expand the function in 
eleven different regions. The formula will be presented in a table format. Each row 
in the table, which represents one region of the expansion, provides the following 
information: 

1. A short-hand label, Ui, for the expansion of u(j],p,q,A) in region 

2. The sign to be used when adding the term to all others. 

3. The range of values that describes the region. 

4. The appropriate expansion parameters. 

Table 14.11 shows how this format would look for the previous case of three variables. 
For a function of four variables, we have the terms listed in Table 14.21 Once again, 
we offer no formal proof of our claim. 

One thing that becomes trickier when working with four variables is knowing which 
terms to keep at any given order. Consider the region <p < g < A, and assume 
we want the expansion accurate to order Oijf / A 2 ). Obviously, first order terms such 
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Label 


Sign 


Region 


Parameters 


Ml 


+ 


77 < (p, g, A) 


77/p, 77/g, 77/A 


u 2 


+ 


(77, g) < (p, A) 


77/p, g/p, 77/A, g/A 


U3 


— 


77 < g < (p, A) 


77/p, 77/g, 77/A, g/p, g/A 


«4 


+ 


(77, p) < (g, A) 


77/g, 77/A, p/g, p/A 


u 5 


+ 


(?7,p,g) < A 


77/A, p/A, g/A 


u 6 


— 


(77, p) < g < A 


77/g, 77/A, p/g,p/A, g/A 


u 7 




77 <p < (g, A) 


77/p, 77/g, 77/A, p/g, p/A 


u 8 




77 < (p, g) < A 


77/p, 77/g, 77/A, p/A, g/A 


u 9 




(77, g) < p < A 


77/p, 77/A, g/p, g/A, p/A 


U10 


+ 


7/<p<g< A 


77/p, 7// g, 77/ A, p/ g, p/A, g/A 




+ 


r/ « ? « p « A 


77/p, 77/g, 77/A, g/p, p/A, g/A 



Table 4.2: Expansion terms for a function of four variables. 



as 77/g, p/A, and 77/A are kept. All products of first order terms that reduce to first 
order terms are also kept (e.g., ij]/p){p/q) = rj/q or (77/p) (p/g) (g/A) = 77/A). 

But what about (77/p) (g/A)? It does not reduce to a first order term, so is it a 
second order term? Is it the same order as (77/g) (p/A)? To answer these questions, 
we present the following rule for determining the order of a product: 

• In whichever region the expansion is being computed, replace <C with <. Mini- 
mize the term subject to these new constraints. If this minimum value is greater 
than or equal to the desired order, the term must be kept. 2 

In this example, the region 77 p q A leads to the constraint 77 < p < q < A. 
A term such as (rj/p)(q/ A) would be minimized by setting p = q. The resulting 
minimum, 77/A, is of the order we want. This implies that any value of this product 
is always greater than or equal to the desired order and must be kept. 

On the other hand, the product (77/g) (p/A) is minimized by setting q = A and 
p = 77. This minimum is 0(r) 2 /A 2 ). However, the product is maximized by setting 
p = q, giving a product of 77/A. So we see that the order of this product is always 
bounded by Oijf /A 2 ) and 0(r)/A). Because its order at any point is less than the 
order we are expanding to, it may be safely dropped. Keep in mind though that the 
product would be necessary for a calculation to 0(r] 2 /A 2 ), which is why we classify 
the product according to its minimum value. 

2 This rule is only needed in two regions: 77 <C p <SC q <SC A and rj <SC q <C p <C A. 
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4.4 Integrating a Four Parameter Function 



Just as we were able to integrate f(r),p,A) to obtain a value accurate to the same 
order as f(r), p, A), we can integrate u(rj, p, q, A) to obtain a function of three variables 
that is equally accurate. Let 

roo 

/(?7,p,A)=/ dqu(r],p,q,A). (4.24) 
J o 

We can now relate the functions fi, fd, and f\ to the integration of the appropriate 
terms of u. fi is determined in the region (rj,p) <C A. Only the terms -u 4 , u 5 , and u Q 
satisfy this requirement. Therefore, 

rOO 

fi(v,P, A ) = dq [u 4 (r],p,q,A) + u 5 (r],p,q,A)-u 6 (r],p,q,A)]. (4.25) 
Jo 

Likewise, 

roo 

f h (r],p,A)= dq [u 1 (rj,p,q,A) + u 2 (rj,p,q,A)-u 3 (rj,p,q,A)], (4.26) 
Jo 

roo 

f d (r),p,A) = / dq [u 7 (ri,p,q,A) + u 8 (r],p,q,A) + u 9 (ri,p,q,A) 
Jo 

- u 10 (r],p,q,A) -u u (ri,p,q,A)]. (4.27) 

As before, the possibility exists that the integration of any one term may be divergent, 
even though the entire quantity is finite. Fortunately, the same subtraction scheme 
defined for the case of three variables can be applied here as well. The result is 
that, for certain classes of functions, some terms in the expansion of u(r],p, q, A) may 
be entirely removed. In particular, the terms u 3 , w 6 , u w , and tin are expanded in 
regions where q is either much larger or much smaller than every other parameter. If 
the function is such that q only occurs in the form of pure powers in the expansion, 
the subtraction scheme demands that these terms be dropped. This leaves us with 
the simplified formulas: 

roo 

fi(v,P,A) = / dq [u 4 (r],p,q, A) + u b (r],p, q, A)] , (4.28) 
Jo 

roo 

f h (rj,p : A) = / dq [u 1 (rj : p,q : A) + u 2 (V:P,q,A)}: (4.29) 
Jo 

roo 

f d (r],p,A) = / dq [u 7 (rj,p,q, A) + w 8 (r?,p, g, A) +u 9 (rj,p,q, A)] . (4.30) 
Jo 

It must be remembered that these equations are not universally applicable, and they 
must be restricted to functions whose expansions for u 3 , uq, u±q, and uu contain 
q only as pure powers. Later, we will see functions that do not allow this type of 
simplification. 
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CHAPTER 5 



THREE-BODY EQUATION EXPANSION 



We are now prepared to apply the technique outlined in the previous chapter to 
the three-body bound-state equation. The expansion will result in two sets of integral 
equations. The first set provides the solution for the leading order behavior. We refer 
to these equations as the LO (Leading Order) equations. The second set provides 
corrections that are smaller than the LO results by a factor of y/B^/ A or \fB~lj A. 
These will be referred to as NLO (Next-to-Leading Order) equations. Because there 
are many steps involved in the derivation of these equation sets, we provide a brief 
overview describing the contents of this chapter. 

First, we begin by setting up some necessary information that will be needed for 
later steps. This includes defining several new quantities and rewriting the bound 
state Eq. (|3.27|) in terms of these quantities. Following this is an expansion of several 
functions that constitute parts of the bound-state equation, and a more detailed 
explanation of how they will be used. 

Next, we use the previously derived results to obtain the equations describing the 
LO behavior. The process is then continued to construct the set of NLO equations. 
Any simplifications or exact solutions to these equations will be shown along the way. 
Finally, we derive the relations needed to obtain values for the three-body coupling, 
g 3 . A summary of all vital equations is given at the end for use as a quick reference. 

5.1 Preparation for Expansion 

To start, we define a few new quantities, including a function that will take the place 
of $(p). This is done to facilitate the expansion later on. Reasons for these new 
definitions will be given as the new quantities are defined. We will also be repeatedly 
applying the expansion techniques in Chapter IH so a basic understanding of these 
methods is necessary to follow the derivation. For readers not concerned with the 
details, a complete summary of all equations is provided in Sec. 15.51 
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5.1.1 Definitions 



Let us begin by denning 



7] 2 = \JB 2 

773 = Jb 3 




(5.1) 
(5.2) 



Using r/2 and 773 in our equations allows us to deal with quantities that have the same 
dimension as the momentum variable. It also identifies their role in the expansion 
as the parameter 77 from the previous chapter. In fact, we will occasionally use 77 to 
generically refer to 772 and 773 in places where either is a valid alternative. 1 
Dimensionless versions of the couplings g 2 and g 3 are defined by 



These dimensionless couplings help to clarify the power counting arguments used 
later. However, for the time being, we choose to replace G 3 with a new coupling 



The reasons for doing so are subtle and not entirely obvious at first glance. The 
quantity 5 is finite and possesses a simple cosine behavior to leading order, whereas 
G 3 will be seen to exhibit divergent behavior at certain points. If we are interested 
only in binding energies, we may solve the integral equation using 5 and need not 
perform the additional integration necessary to compute $1 and hence G 3 . These 
statements will be justified later. 

Another dimensionless quantity, and perhaps the most important, is the replace- 
ment of the pseudo-wavefunction with 



This shares some of the benefits of all the previous definitions. Like 7] 2 and 773, its role 
in the expansion is easily seen to be the same as the function / used in the previous 
chapter. Like G 2 and G 3 , the dimensionless nature of / will simplify future power 
counting. In addition, this function tends to be of 0(1) throughout the entire range 
of p and is less prone to large numerical fluctuations than $ when B 2 ~ B 3 . 3 

1 For example, saying that something is 0(r]/A) can imply that it is (9(r?2/A) or (9(r?3/A). 

2 Although S is not technically a coupling, we will often refer to it as such because it plays a role 
very similar to G3. 

3 We cannot prove these statements a priori and must verify them numerically. 



G 2 = Ag 2 , 
G 3 = A 4 g 3 . 



(5.3) 
(5.4) 



f(rj 2 ,r] 3 ,p,A) = {p 2 + r 1 l-r 1 l)§(p). 



(5.6) 
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With these new definitions, the three-body bound-state equation (|3.27|) now takes 
the form 



5.1.2 Conventions 

The process of expanding Eq. (|5.7|) is done in a series of smaller steps. Because the S 
term is simply added to the integral term, each can be expanded individually and then 
added together at the end. In addition, these terms are composed of other quantities 
like Di(p), G2, etc. which have their own expansions in terms of 77/A. The goal is to 
split the equation for / into separate equations for f], fd, and fh- At the same time, 
we want to distinguish between the LO and NLO portions of these equations. In 
order to categorize the contributions from every function in the three-body equation, 
we develop a labeling convention. 

Just as the function / can be divided into its low-, mid-, and high-momentum 
parts, so too can any other function of 77, p, and A. The same subscripts will be 
used to denote these terms. For example, the approximation to D(p) in the region 
p ~ i) C A will be written as Di(p). This low-momentum version may also be 
expanded in powers of 77/ A and p/A when attempting to uniformly expand D(p) to 
a given order in 77/A. To distinguish between orders, we will use a subscript n to 
represent the 0((rj/A) n ) portion. Thus, Du(p) denotes the part of Di(p) having the 
same magnitude as 77/A and p/A. Likewise, the function D\ho{p) is the 0(1) part of 
Di(p) in the region 77 <C p ~ A. 

The couplings are labeled in a similar manner. Since they contain no p dependence, 
the classification of I, d, or h has no meaning. Yet they may be expanded in powers of 
77/A, so the numerical subscript is useful. For the coupling 5, the order n contribution 
will be identified as 8 n , which should not be confused with the the phase shift 5o(p). 
For a coupling like G2 which already possesses a subscript, we shall write it as Gi, n - 
This comma delimiter takes on the role of '/', 'cf, or '/V, acting as a placeholder 
between the two numbers. This helps to avoid confusion in cases where the subscript 
may already consist of a double-digit number. 

5.1.3 Expansion of Components 

In order to achieve results accurate to 0(r)/A), the following quantities must also 
be expanded to the same order: 5, G 2 , D \ —t}\ — |p 2 ), and Di(p). The value of 5 is 




- 5 



(p 2 + r)l - r)l)D x {p) 



(5.7) 
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unknown prior to solving the equations, so we shall simply refer to the 0(1) part as 8$ 
and the 0(rj/A) part as S\. Ultimately we must show that Sq is indeed 0(1), which is 
equivalent to showing that the 3-body contact interaction is required to renormalize 
the bound-state equation. 

The terms in the expansion of g 2 have already been computed in Eq. (|4.22|) . 
making it easy to write the expansion for G 2 . 

G 2t o = „, 2 ; o; — ~T7i (5-8) 



3/il + 8h 2 + 16 
2048\^tt 2 



° 2 ' 1 ~ (shi + s^l ie) 2 (a) • (5 ' 9) 

Here we have assumed that h 2 is of 0(1), but the above equations will continue to 
hold provided that h 2 does not become large enough to negate the power counting 
arguments. For r e < 0, Eq. ()2.34|) shows that h 2 will diverge like \/A as A — > oo. 
In this case, we would be required to limit the value of A so that h 2 never becomes 
arbitrarily large. 

Next we must expand 

This will be our first application of the methods in Sees. 14.31 and 14.41 We make D fit 
the form of / dq u(r],p, q, A) by defining 

q 2 (\v 2 + U 2 (q) 2 
u( V2 ,r) 3 ,p,q,A)= v J (5.11) 

2tt 2 (2q 2 + r]l) 2q 2 + § p 2 + rj 2 



The function u must now be expanded in each of the eleven different regions up to 
and including first order corrections: 



U2 = ^WTl' (5 ' 13) 
1 

47T 2 



«3 = 73. ( 5 - 14 ) 



8ttV V2 J 
?M 

2vr 2 (2g 2 + % 2 ) (V + \p 2 + % 2 ) 



q 2 (fp 2 + % 2 - V 2 ) 
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8ttV 
16n 2 q 2 



u 7 = J^U 2 (q) 2 , (5-18) 



«8 = -7 ^ = n x , (5-19) 



8tt 2 (2g 2 + §p 2 ) ' 

M9 = ^wra ' (5 - 20) 

-io = Y6ttV' (5 - 21) 
un = JL (5.22) 

The approximations to D in each of the momentum ranges are given by the following 
relations: 

roo 

Di = dq [u 4 + u 5 - u 6 ] , (5.23) 
Jo 

roo 

D d = / dq [u 7 + u 8 + u 9 - u w - u n ] , (5.24) 
Jo 

roo 

D h = / dq [u-y + u 2 - u 3 ] . (5.25) 
Jo 



Now we must differentiate the LO behavior from the NLO behavior for each of the 
above integrands. To do this, we integrate each term and look at the final result 
to determine the order at which it contributes. This is much simpler than trying 
to determine a priori which integrands are LO or NLO and then performing the 
integrations, although we will see some cases where this is necessary. The results are 
shown below. 



n 1 foo q 2 (IP 2 + Vl ~ Vl) 



2tt 2 Jo ( 2 g 2 + rf 2 ) (V + |p 2 + % 2 

\l\p 2 + ri -r]2 

8V2tc 
\p 2 + V 2 ~ V2 f 00 dq 



(5.26) 



Ai(?72,?73,P,A) = 5-3 / — U 2 (q) -1 

oir z Jo q l 

fh 2 2 + 8h 2 -16\ (lp 2 + r]l-r]l 



64v / 2W 2 J \ A 

2 i-OO 



(5.27) 



16tt 



(5.28) 



41 



D dl (r) 2 ,p,A) 



dq 



V2 



\Z 2 q 2 ^ q)2 ^ + \2tt 2 (2g 2 + V 2 ) 
3(h 2 2 + 8h 2 - 16) p 2 



An 2 



+ 

dq 



128v / 2vr 3 / 2 A 

UM 2 



V 



2q 2 + \ V 2 

^^K-3V 2 ) 2 (l-Erf(g)^ 



256tt 2 A 4 

+ 3V2nh 2 Ap ((8 + h 2 ) A 2 - 3/i 2 ;p 2 ) 
g 2 1 



cxp 



(5.29) 

2A 2 ) 
(5.30) 



/ 



2tt 2 (2g 2 + rf 2 ) An 2 



These same steps are used to compute the expansion of Di(p). In 
function u takes the form 

q 2 f 1 , U 2 {q+\p)U 3 {q)U 3 {p)U 3 {p + q\ 



Q f 

u(r} 3 ,p,q, A) = — J dz 



2p 2 + 2q 2 + 2pqz + rg 
and yields the following expansions in each of the eleven regions: 

,2 



(5.31) 
this case, the 

(5.32) 



Ml = 

U 2 = 

U3 = 

U4 = 

U 5 = 

u 6 = 

u 7 = 

u 8 = 

Ug = 



Q 



Atx 2 J-l 

J2 



1 dz \p)UMUs(p)Us(p + q) 



2p 2 + 2g 2 + 2pqz 



q 



U 3 (p) 2 U 2 (p/2), 
U 3 (p) 2 U 2 (p/2), 



Aix 2 p' 

q 2 
Ai\ 2 p 2 

^- 2 U 3 (q) 2 U 2 (q), 

q 

8n 2 p 
1 



In 



rjl + 2p 2 + 2g 2 + 2pq 
j) 3 + 2p 2 + 2g 2 — 2pq 



Air 2 ' 
1 

A^ 2 

q 

8ix 2 p 

q 2 
Aii 2 p 2 ' 



U 3 (q) 2 U 2 (q), 



In 



' p 2 + q 2 + pq s 
p 2 + q 2 -pq, 



(5.33) 
(5.34) 
(5.35) 
(5.36) 
(5.37) 

(5.38) 
(5.39) 
(5.40) 

(5.41) 
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An 2 ' 



4tt 2 p 2 



(5.42) 
(5.43) 



Once again, the LO and NLO terms are identified after integration. It is easily verified 
that the appropriate expansions are 



D m {A) 



D ld0 (A) 



D ldl {p) 



dqU 3 (q) 2 U 2 (q) 



47T 2 JO 

A (6 {h\ + Ah 3 + 12) + h 2 (5h 2 + 12h 3 + 12)) 



In 2 jo 



dq 



576^3 ?r 3 / 2 
q ( yj + 2p 2 + 2q 2 + 2pq 

v yvi + ^p 2 + 2q 2 — 2pq, 



71 



-^/2 V 2 + 3p 2 , 



47T 2 JO 



dqU 3 (q) 2 U 2 (q) 
A (6 (h 2 + Ah 3 + 12) + h 2 (5/i§ + 12h 3 + 12)) 



47T 2 JO 



dq 



576^^3/2 

q . ( ' p 2 + q 2 +pq s 
— in 



(5.44) 



(5.45) 



(5.46) 



2p 



p z + q z — pq 



1 q^ _ 

p 2 p 2 



71 



D lha (p, A) 
D lhl (p) 



2 



P, 



An 2 Jo 
0. 



dq q 2 



1 d U2(g + \p)UMUz{p)Mp + q) 

-i Z 2p 2 + 2q 2 + 2pqz 



(5.47) 

(5.48) 
(5.49) 



Notice that is in some ways easier to analyze and in other ways harder. Because 
u 2 = u 3 , the expression D lh = J dq [u\ + u 2 — u 3 ] reduces to an integration of only 
one term. Unfortunately, this integral has no simple analytic solution. Dimensional 
analysis dictates that the final result must have dimensions of momentum, and since 
the only quantities involved in the integral are p and A, the answer must be O(A) 
or 0(p) = O(A). This is the same order as both Du and -D^o, proving that this 
integral is a LO term. This also means that there are no integrands that contribute 
to NLO, which is why -Dihi = 0. 



5.1.4 Expansion of Main Integrand 

We are now in a position to expand the main integrand of Eq. ()5.7|) . If we ignore the 
three-body interaction for a moment, the bound-state equation can be written as 

POO 

f(V2,V3,P,A)= dqQ(r] 2 ,r] 3 ,p,q,A) f(r] 2 ,r) 3 ,q,A), (5.50) 
Jo 
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where 



nf A x <f fp 2 + vi-vi 



2n 2 D (- % 2 - \p 2 ) W + Vl - Vl 



x /' (5.51) 

By equating u(rj 2 ,r] 3 ,p,q, A) = Q(r] 2 , i] 3 ,p, q, A) /(r^, 773, <?, A), the bound-state equa- 
tion becomes /(t7 2 , 77 3 ,£>, A) = J dq u(r] 2 , r] 3 , p, q, A). This is just the form needed to 
expand the integral. 

The function Q can be expanded in the same 11 regions as u. On the other 
hand, the function / depends only upon three scales: 77, p (or q), and A. It is 
expanded in only three regions to produce fi, fa, and fh- The relationship between 
the expansions of u, Q, and / is easily determined by paying attention to the region 
of the expansion. For instance, u\ results from an expansion where 77 <C (p,q,A), as 
does Qi. In this same region, f(q) would be expanded to fh(o)- Thus ui(r),p, q, A) = 
Qi(i],p,q,A)f h (r],q,A). By th is sam e reaso ning, u 6 = Q 6 f d (q), u 9 = Q 9 fi(q), etc. 

With this in mind, Eqs. (|4.25|) . (j4.26|) . and ()4.27j) can be used to obtain the 
following relations: 

flip) = / dq[Qtf h (q) + Q 5 Mq)-Q 6 f d (q)), (5.52) 
Jo 

fd(p) = / dq [Q 7 f h (q) + Q 8 fd(q) + Q9fi(q)-Qiofd(q)-Qiifd(q)], (5.53) 
Jo 

HP) = / dq[Q 1 f h (q) + Q 2 f l (q)-Q 3 f d (q)]. (5.54) 
Jo 

To avoid notational clutter, we have dropped the explicit dependence of the / func- 
tions on i]2, 773, and A. We will continue to do so unless confusion would result or 
such dependence needs to be emphasized. 

The task of expanding the integral equations now becomes the task of expanding 
Q. Because Q contains the function D (—rfe — |p 2 ) , each Qi will be written in terms 
of the expanded D functions ()5.26|) through ()5.31|) . The eleven expansion terms are 
shown below. 



Q, = , Jl - ^pl) r it U > (f + (f + & (5.55) 



27T 2 ^o(p,A)V D h0 (p,A)JJ-i ' 2p 2 + 2q 2 + 2pqz 
Q2 = 9 T( 2+ 2 ^ 2s n 1 A , f 1 - n hl ( {r]2 l ) U2(p)U,(p/2) (5.56) 

Q 3 = - 2n \ tt fl- ° hl ( {V2 l )u 2 (p)U2(p/2) (5.57) 
2n 2 D h0 (p,A) \ D h0 (p,A)J 

(p 2 + r)j - rj) ( D n (ri2,ri3,p,A) \ TT ^ TT 

Qi = 2 2 n 7 r 1 7^-7 x- C/2(g )U 2 (q 2) 5.58 

27v 2 q 2 D l0 (ri 2 ,r] 3 ,p) \ D l0 (r] 2 ,T] 3 ,p) J 
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i L _ Dn(v2,V3,p, A ) \ q (p 2 + vi~ gf) 

4:Tr 2 Dio{r}2,rj 3 ,p) \ D l0 (r) 2 ,r) 3 ,p) J p (q 2 + rfc - T] 2 ) 
'77J + 2p 2 + 2g 2 + 2pg 
7$ + 2p 2 + 2g 2 - 2pq 



x In ' /3 2 - *JL - ^ 9 - ^ (5.59) 



2n 2 q 2 D l0 (r]2,V3,p) \ Aofe,%,p) / 

2n 2 q 2 D d0 {p) \ D d0 (p) J 

n p (1 D di(V2,P, &) \ 1 ( p 2 + q 2 +pq \ , „ Q s 

g§ - i^w^i 1 - ^ o(P ) rU+^-^J (5 - 6 ) 

^ 9 27r 2 (g 2 + 77|-77|)^ (p) V A»(p) ) ^ 

Qw = 2^DM\ l (5 - 64) 

On - (5-65) 
2Tc 2 D d0 (p) \ D d0 {p) J 

Of course, the three-body interaction has been neglected in Eqs. ()5.52j) to (|5.54jl . and 
it must likewise be expanded and included. 



5.2 Leading Order Equations 

Most of the necessary work has been completed, allowing us to begin expanding the 
integral equations and separating the LO and NLO behavior. The leading order 
is handled in this section, with next-to-leading order corrections being treated in 
Sec. O 

Discerning the order at which certain terms contribute now becomes more difficult. 
Unlike previous expansions, we cannot integrate the terms to see the behavior of the 
final answer. The reason is that some of the integrals involve /, which is the unknown 
function we are trying to solve for. Therefore, we must resort to dimensional analysis 
and other methods that may require knowledge of the final solution. We provide a 
detailed explanation of each technique the first time it is used, but the details are 
omitted for later uses in the interest of brevity. 

5.2.1 Low-Momentum Region 

Consider the expansion terms comprising fi(p): 

U 2 (q)U 2 {q/2)f h (q) 



MP) 



p 2 + vi 



■nl 



2n 2 Dio(r] 2 ,ri 3 ,p) 



Dn{ri2,ri3,p,A) 

Dlo(V2,V3,P) 
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+ 



2p(q 2 +rg- r$) 




(So + 5i) 



(P 2 + Vl - vl)(D m (A) + D m ( m ,p)) 



(5.66) 



A 2 (G 2 ,o + G 2A ) (D l0 ( V2 , m ,p) + D n (rj 2 ,rj 3 ,p,A))' 



For the leading order equation, we need to identify only the terms that are 0(1). 
Obviously, the ratio Dh/Diq can be dropped from the factor in front of the integral 
because it only provides arj/A correction compared to 1. The remaining part of this 
factor behaves as 0(r)), and the integral itself must have the dimensions of inverse 
momentum to make the product dimensionless. 

Next, we analyze each integrand to determine its leading order behavior. The 
first integrand contains the functions U 2 and fh- For a leading order analysis, we only 
need to consider fho(q)- The dimensionful parameters in this function are q and A. 4 
The same holds true for U 2 . This means that after integrating over q, only A will be 
left to provide the needed dimension to the result. Consequently, the integral must 
behave as 0(A~ l ). With the factor in front of the integral, we are left with a term of 
0(r)/A). Thus, we can drop this term for now and use it later when we analyze the 
NLO behavior. 

Before continuing, notice that this integrand possesses a divergence as q — > 0. 
In this limit, U 2 (q) approaches 1, and fho(q) approaches fdo(q)- We will see later 
that fdo is proportional to a cosine, making the integral diverge like 1/q. This is 
not a problem because the fdo term from the third integrand cancels this divergence. 
Together they yield a finite result. This allows us to classify the third integrand as a 
0(r]/A) contribution and drop it from immediate consideration. 

The second integrand involves a term containing r] 2 , i] 3 , and p along with the 
function fi(q). Again, only the Oil) part of namely fio, needs to be considered. 
As we will soon see, it has no A dependence. Therefore, the integral only involves 
quantities of order r\ and must behave as 1/rj. Combining this with the factor in front 
of the integral, we get the only 0(1) contribution from the integration. 

The last part to investigate is the three-body interaction term. If we take only 
the leading order parts of each piece, we find that the entire term behaves like 



This can only make an addition to the NLO behavior. 

Overall, we find that only the second integrand is used in the equation for fi . 
The final equation is 




(5.67) 



fio(r)2,V3,P) 



4n 2 Di (r] 2 ,ri 3 ,p) Jo 




q (p 2 + rj) 
p (q 2 + r] 2 - 7%) 




(5.68) 



This statement will be verified later when we examine the equation for fho- 
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From this equation, we see that our previous assumption that fi contains no A 
dependence is indeed correct. 



5.2.2 Mid-Momentum Region 

In the mid-momentum region, the integral equation is 



1 (, D dl (ri2,p,A)\ f°° \p 2 



2ix 2 D d0 {p) \ D d0 {p) J Jo [q 2 

p ( p 2 + q 2 + pq \ q 2 p 2 

+ ^r ln 2 ; 2 M?) + 7-5-; — 2 — - - M?) 

2g + q 2 — pq J (q 2 + ^3 — rji) q 2 

, r 1 s) (P 2 + Vl ~ V 2 2 )(Dido(A) + Didijp)) ( , m , 
{ 0+ 1J A 2 (G2,o + G 2 ,i)(^o(p) + ^i(r/2,P,A))' l ' J 

The same reasoning used for / w can be applied here as well. The D dl /D d0 ratio is 
dropped from the integral prefactor, resulting in a 1/p term. The first integrand has 
a q —>■ divergence that is canceled by the fourth integrand. Together they yield a 
finite integral that dimensional analysis shows to be 0(p 2 /A). Combined with the 
prefactor, the end result is a 0(p/A) term that does not contribute to leading order. 

The third integrand is divergent as q — > 00, but the fifth integrand cancels this 
divergence. The result is a finite term that behaves as rj/p, also a NLO correction. 
Only the second integrand contributes to leading order. The function f d (q) is replaced 
by f d o{q), and the integral of this term is of order p. With the prefactor, the term is 
Oil). 

The leading behavior of the three-body interaction is 

OWJ1 ^L-o(^, (5.70) 



A 2 0(l)0(p) \A 

which cannot contribute to leading order. Thus, only one integrand adds to the 
leading order equation for f d , giving us 

t l \ I r dq ( p 2 + g 2 + pq \ 

fdoip) = -: 2n pr / —In fdo{q). (5.71) 

4n 2 D d0 (p) Jo q \p 2 + q 2 - pq J 

This equation happens to have an analytic solution [T2] . 

f d0 (p) = A cos (s In (I) + 0) , (5.72) 

where A is an arbitrary amplitude, 6 is a phase determined by boundary conditions, 
and s is the real, positive solution to the equation 
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The value of s is approximately 1.006237825102. 

The A dependence in this solution may be a bit misleading since Eq. (|5.71|) does 
not contain it. The cutoff is included so that the argument of the logarithm is 
dimensionless. Making this choice is a matter of preference, and it could just as 
easily have been 773. 



5.2.3 High-Momentum Region 

The expansion terms for the high-momentum region take the form 



hip) = 2n 7 — v\ / dq 

2iT 2 D h0 (p,A) 







r 



P Li 2p* + 2q 2 + 2pqz Mq) 



+ U 2 (p)U 2 (p/2) ^ —^—^ /,(g) - f d (q)^ 

_ , r rs + -nl - ri)D lhQ {jp, K) 

1 + A 2 (G 2i0 + G 2jl ) (D h0 (p, A) + D hl ( V2 )) ■ 1 • 1 

For leading order behavior, we drop the ratio D^i/D^ from the integral's prefactor. 
The remaining part of the prefactor behaves as 1/A. The only parameters in the first 
integrand with the same dimension as momentum are p and A, both of which are of 
the same order. Dimensional analysis requires that this integral be 0(A), making its 
product with the prefactor 0(1). Thus, it contributes to the leading order equation. 

The third integrand cancels the q — > 00 divergence in the second, resulting in a 
finite term of order r//A as required by dimensional analysis. These terms are NLO 
contributions that we will need later. 

Last of all, we examine the leading order behavior of the three-body interaction. 
We find it to be 

p 2 D lh0 (p,A) p 2 Q(A) 

- *° a^wwy) ~ vowm ~ (5 - T5) 

making it necessary to include it in our equation for fho- This is the only leading 
order equation that contains a three-body interaction term. The exact form of this 
term is 

3/^ + 8/^ + 16 ( D lh0 (p,A) \ (p 2 \ 
" 60 1 2 8v /2V/ 2 { A^AyJ [T 2 ) ' (5 - 76) 

giving us the last piece we need to construct the complete integral equation for fho- 

Ao(p,a) = Tmfi rr / <k dz — o 9 , o % , o L fm{q,A) 



2ir 2 D h0 (p, A) Jo J-i 2p 2 + 2g 2 + 2pqz 

3/i 2 + 8/i 2 + 16 ( D lh0 (p,AY ,, , 
5 ° 12872^ lA^AT 11 - 1 - (5 - 77) 
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Notice that this equation, and hence the solution /^ , only contains dependence on 
p and A. This fact was used earlier in our derivation of the integral equation for fi , 
and now it has been justified. 



5.3 Next-To-Leading Order Equations 

The leading order equations for all three pieces of / have been derived. Now we 
must go back and gather the NLO terms that were neglected in our first pass, as 
well as determine any sub-leading behavior to the leading order terms that may also 
contribute. However, before doing so, we derive two relations that will be needed for 
our first order corrections. 



5.3.1 Leading Order Identities 



As p becomes much greater than r], the function fmip) must begin to equal fdo(p)- 
If we consider the leading order behavior of Eq. (|5.68|) when p ^> rj, we obtain the 
relation 



fdoip) 



4irW d0 (p) Jo 



dq 



qp 



p 2 + q 2 + pq 



q +T}3 — r\2 \p + q — pq 



In 



fio(V2,V3,q)- (5.78) 



Since f d o(p) satisfies Eq. ()5.7H) . the above relation can also be rewritten as 



Au 2 D d0 (p) Jo 



dq In 



' p 2 + q 2 + pq 



p2 _|_ g2 



PQ, 



qp 



P 



q 2 + rjl 



fio(V2,V3,q) - -fdo(q) 



q 



-- 0. 

(5.79) 

For values of q ^> rj, the integrand will vanish. Therefore, the previous relation must 
also hold when we restrict the q integration to values of the same order as rj. Because 
p is much greater than rj, and hence q, the logarithm can be approximated as 2q/p. 
The relation now becomes 



47r 2 D d0 (p) 
or more simply 



Jo p 



qp 



P 



dq 



q 2 + m 



vi 



fio(V2,V3,q) - -fdo(q) 



0. 



q + m - m 



/zo (772,%,?) - /do (q) 



0. 



(5.80) 



(5.81) 



Just as fio(p) equals fdo(p) for p ^> r], so must fhoip) for p A. The leading 
order behavior of Eq. (J5.77|) for this range of p is 



fdoip) 



P 



27rW d0 (p) Jo 



dq 



-1 2p 2 + 2q 2 + 2pqz J UW ' ' 



Si 



3h 2 + 8h 2 + 16 ( D ldD (A) \ (p^ 
128^^/2 I D d0 (p) ) \A 2 / 



(5.82) 
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Using Eq. ()5.71|) once again allows us to rewrite the previous equation in the form 

= TTTTT^ r-^f 4x4 j ^ )^(g)W 2 )Ao(g,A)-/ rfo (g)] 
4ir 2 D d0 (p) Jo q \p 2 + q 2 - pq J 

~ 6 ° 128^^ {-JJMJ \A 2 ) ' ( } 

The integral vanishes for values of q much less than A, so it may safely be restricted 
to values greater than or on the order of A. This implies that p <C q, permitting us to 
approximate the logarithm by 2p/q. Using this approximation gives us the integral 

r <lq [U 2 (q)U 2 (q/2)f h0 (q,A) - f d0 (q)) . (5.84) 



2TT 2 D do (p)Jo q 2 

After dividing out the integral's p-dependent prefactor, we are left with the relation 

= r^[U 2 (q)U 2 (q/2)f h0 (q,A)-f d0 (q)] 
Jo q z 



So g4 ^ A2 + 8^2 + 16J Dido(A). (5.85) 

Equations (|5.81|) and f|5.85|) are the key relations, and in a moment we shall see why 
they are necessary. 

5.3.2 Low-Momentum Region 

As shown before, the expansion for f] is 

P 2 + vl~V2 D n (rj 2 ,rj 3 ,p,A)\ [°° dq 



MP) 



2ir 2 D l0 (r) 2 ,r) 3 ,p) \ D m (rj 2 ,rj 3 , P) J Jo Q 



,2 



U 2 (q)U 2 {q/2)f h (q) 



, (f 1 ( v! + 2p 2 + 2q 2 + 2pq \ 

Mq 2 + vi - vi) Ui + V + V - 2p? J Jl{q> Jd{q \ 

k , n (P 2 + Vl - V 2 2 )(Duo(A) + Dm(m,P)) 

l0+ l) AZ(G 2 , + G 2 , 1 )(D l0 ( V2 , V3 ,p)+D ll (r l2 ,r l3 ,p,A)y { ' ) 

The leading order behavior has been extracted to get Eq. ()5.68|) . and now we wish to 
examine the terms of order 77/A to obtain the NLO equation. 

As noted is Sec. 15.2.11 the leading behavior of the first and third integrands is 
0(r]/A). Together, these terms provide an addition to the equation for first order 
corrections: 

,r + r (Ul [U 2 (q)U 2 (q/2) f h0 (q) - f d (q)} (5.87) 



2ir 2 D l0 (ri 2} ri 3 ,p) Jo q 
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Using (|5.85p . this may be simplified to 



p 2 + vi 



V2 



2n 2 Di (r] 2 ,r] 3: p) 



^2 „2' 



(p 2 + vi-V 2 2 )D ld0 (V 
A 2 G 2fi Dio(m,V3,p) 



(51 



Any other pieces of the first and third integrands are necessarily of higher order and 
may be ignored. 

The leading order behavior of the second integrand was determined to be 0(1), 
which means there may be sub-leading terms that are of order 0(r]/A). In fact, there 
are two such terms. One results from dropping the Dh/Diq ratio from the integral's 
prefactor and using fu in the integrand. The other retains the prefactor ratio and 
uses fio in the integrand. Their combined contribution is 



p 2 + r] 2 



V2 



4iT 2 pDi (ri2 } ri 3 ,p) 
Dn(m,V3,P, A) 



<r + vi 



In 



dq 



q 



In 



(r] 2 + 2p 2 + 2q 2 + 2pq 
\rii + 2p 2 + 2q 2 - 2pq 
T)l + 2p 2 + 2q 2 + 2pq\ 



Ao {rj2, Vs, P) Jo q 2 + r/f - r)\ \T% + 2p 2 + 2q 2 - 2pq 




(5.89) 



which may also be written as 



4ir 2 D l0 (ri2,r] 3 ,p) Jo 



dq 



q(p 2 + 7]l 



p{q 2 +vi~V2 




ji + 2p 2 + 2q 2 - 2pq 
_ Dii(V2,V3iPi A) 

Dl (V2,V3,P) 



fioip) 



(5.90) 



by using Eq. ()5.68|) . 

The final source of possible corrections is the three-body interaction term. Our 
previous analysis of fa shows that the leading behavior of this term is 0(j]/A). There- 
fore, it contributes a term to the NLO equation that looks like 



A 2 G 2 ,oD lo (r)2,V3,p) 
Combining the results of (jSHHj) , (JSHOI), and (l5~TJTl) yields 



(5.91) 



A-K 2 D l0 (r] 2 ,7] 3 ,p) Jo 
Dn(v2,V3,P, A) 



dq 



Mo) 



Dio{r]2,r} 3 ,p) 



fio(p) + So 



~ q(p 2 + r£ - rj) ( rj + 2p 2 + 2q 2 + 2pq 
p(q 2 + 7] 3 — 7] 2 ) \r) 3 + 2p 2 + 2q 2 — 2pq / 
\p 2 + % 2 - r, 2 )D im {A) x (p 2 + V 2 - ri)D m {A) 



A 2 G 2fi D m {r] 2 ,rj 3 ,p) 



-So 



A 2 G 2 ,oDio{r]2,ri3,p) 



which may be simplified by recognizing that Dido{A) = Duo(A). The final equation 
for the first order correction in the low-momentum region becomes 



fii{r} 2 ,r} 3 ,p,A) 



1 



4ir 2 D l0 (ri 2} ri 3 ,p) Jo 



dq 



q(p 2 



p(q 2 



vi 



vi) 
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f rj + 2p 2 + 2q 2 + 2pq . 
Ai(%,?73,P, A) 



Dio(v2,m,p) 



fio(V2,V3,P)- (5-92) 



5.3.3 Mid-Momentum Region 



The first order correction terms for /^i are obtained in a manner very similar to 
/tt.. The leading order analysis of Eq. (|5.69|) showed that the first, third, fourth, and 
fifth integrands have an 0(r)/A) behavior for their leading terms. This provides a 
contribution of 



dq 



p 2 p 2 q 2 
-^U 2 (q)U 2 (q/2) f h0 (q) - -^f d0 (q) + -5-; — 2 5/10(0) ~ /do(?) 



2n 2 D d0 (p) Jo 

The identities (|5.81|) and f)5.85j) can be used to reduce these terms to 



(5.93) 



p 2 D ld0 (A) 
A 2 G 2 ,oD d o{p) 

The three-body interaction also provides a term of the same order, 

p 2 D ld0 (A) 
A 2 G 2 ,oD d Q{p) 

which exactly cancels the previous integral terms. The only term left to provide 
corrections is the second integrand involving f d (q). There are two obvious corrections: 
one that results from dropping the D d i/D d0 ratio in the prefactor and using f d \(q) 
in the integral, and the other that results from keeping the ratio in the prefactor 
and using f d o(q) in the integral. However, there is a third contribution that involves 
keeping the prefactor ratio and using f d i(q) in the integral. For the case of fn, no 
such term is encountered because the product would be Oirf / A 2 ). Yet, in this case, 
it is entirely possible that two first order corrections can have a product that is also 
a first order correction. 5 

The final integral equation is now written as 



f i x V L D dl (7] 2 ,p,A)\ r°° dq p 2 + q 2 +pq\ (f , x . , M v 

fdi{P) = Y~2~n~T~\ 1 n~~T\ — / — ln 2 7 2 (/*>W + w) ' 

ATT 2 D d0 (p) \ D d0 (p) J Jo q \p 2 + q 2 -pqj 

(5.96) 

and we must remember to keep only the terms that are the same order as f d \. 
Like the leading order function, f d o, this equation has an analytic solution: 



f d i(r)2,p,A) = An— cos ^s ln (j^j +8^ 



5 



The simplest example is (f]/p)(j>/ A) = (77/A). 
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+ A hl ^cos (s hi (j^ 

A V2 ( , fP 

°X COS I VA 



~0 W 

+ ddlO 



+ A dll ^ln Ij^J cos ^s Q \n ( ^ ) i 



|) cos (so In (| 

The details of the derivation leading to this form are given in Appendix 1X1 



(5.97) 



5.3.4 High-Momentum Region 

During the analysis of Eq. (J5.74j) to determine fho, we discovered that the second and 
third integrands have a leading behavior that can only contribute to the first order 
corrections. Upon closer inspection, we see that the sum of these integrands fits the 
form of Eq. (|5.81|) and must therefore equal zero. 

The first integrand provides two additions to the NLO result, 



P 



and 



2nW h0 (p,A) Jo 
p 2 D hl (r] 2 ) 



dq 



l d U 2 {q+\p)U 2 (p+\q) 
-l 2p 2 + 2q 2 + 2pqz 



,1 lh (q+ W) lh (P+ \q) 
2^D m ( P ,kyh a<1 Li 2p 2 + 2q* + 2pqz Jm{qh 

Using the integral equation for fho(p), equation ()5.99|) can be simplified: 



Dhi(m) 
D h0 (p,A) 



ho{p) + 5 



p 2 D lh0 (p,A) 



(5.91 



(5.99) 



(5.100) 



A 2 G 2 , D M (p,A) 

Finally, we must consider the three-body interaction term. The leading behavior 
of this term is 0(1), but there are three sub-leading terms that must be added to our 
first order corrections. These terms are 

p 2 D lh0 (p, A) 



-ft 



A 2 G 2 ,o D h0 (p,A)' 



and 



& 



So 



p 2 D lh0 (p,A) G 2 ,i 
A 2 G 2 2fi D h0 (p,Ay 

p 2 D lh0 (p,A)D hl (r] 2 ) 



(5.101) 
(5.102) 

(5.103) 



A 2 G 2fi D hQ ( Pl A) 2 ' 

Notice that Eq. (J5.103j) cancels the second part of Eq. (J5.100)) . The final integral 
equation is 

-I . U 2 (q+\p)U 2 (p+\q) 



fhl(V2,V3,P,&) 



P 



2ir 2 D h0 (p,A) Jo 
D h i(r] 2 



dq 



dz ■ 



2p 2 + 2q 2 + 2pqz 



D h0 (p, A) 



ho(p, A) - 5 



p 2 D lh0 (p,A) 
A 2 G 2fi D h0 (p,A) 



Si G 2 ,i 

To 



G 



2.0 



/w(??2,^3,g, A) 
(5.104) 
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5.4 Obtaining Values for G 3 



With the integral equations derived so far, it is possible to compute solutions for all of 
the / functions as well as the value of 5. The next step is to determine the necessary 
equations so that the value of 5 can be used to compute the value of G 3 . The basic 
equations needed are: 

6 = G 3 $ 1 /A, 



3G 9 



2tt 2 A 

By defining the dimensionless quantity 

3 



dp p 2 -Di(p) 



G 3 D 2 ^ 

A 4 ' 



1 



dp p 2 -Di(p) $(p) 



2vr 2 A 2 

the equation for $i may be written as 

$i = G 2 A1-G 3 D 2 $ 1 /A\ 

or equivalently 

6A -CA1 D * 6 

The goal is to expand both sides of this equation in powers of 77/ A and then equate 
terms of the same order on both sides. The result will be a relationship between 5q, 8\ 
and G 30 , G 3 1. From there, the value of g 3 is easily reconstructed. 



(5.105) 
(5.106) 

(5.107) 

(5.108) 
(5.109) 



5.4.1 Matching Conditions 

We know the expansion for every quantity in Eq. ()5.109j) except for two: I and D 2 . 
We will expand X in a moment, but for now we simply write the terms of its expansion 
as Iq and T\. The remaining quantity is defined by 

D2 = jl r> p r fdq f 1 dz uj^wmM+P? (5 110) 

87r 4 io y y Jo H H J-i 2p 2 + 2q 2 + 2pqz + rfc V ; 

The order (r/ 3 /A)° term is obtained simply by setting 773 equal to zero. This results 
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D2 ' = 8^io Mp io 9d9 l/ z ^WT^T' (5 ' ni) 

which by dimensional analysis must be proportional to A 4 . Hence, D 2 ^/A 4 is a 
constant. The next term in the expansion of D 2 behaves like (773/A) 2 , so we find that 
£2,1 = 0. 

If we multiply both sides of Eq. ()5.109|) by G 3 /A, we obtain 

5 = G 3 (G 2 l-^5y (5.112) 
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By expanding both sides, 

So + 5i = (G 3)0 + G 3j i 
and matching corresponding orders, 
= G 3)0 ^G 2) o X <5 C 

Si = ^3,0 f C?2,0 ^1 + C?2,l ^0 

we obtain the relations 

G3.0 



(G 2 ,o + G 2il )(Jo+Xi; 
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2.0 



A 4 



(So + Si 



D 



2.0 



A 4 



Si + G 3>1 G 2 , 1( 



D 



2.0 



A 4 



Sn 



G 2fi l -D 2fi S /A*- 



G 3 A 



Si 



1 



Sn 



On 



^2,0 

A 4 



3,0- 



(5.113) 

(5.114) 
(5.115) 

(5.116) 
(5.117) 



5.4.2 Expansion of X 

To expand X, we utilize the expansion technique of Sec. 14.21 The definition 

3p 2 



2tt 2 A 2 G» 2 + rft -r)l 



-Di(p)f(p) 



(5.118) 



lets us write X = J dp[ui + Uh — u d ]- The function u must be expanded to 0(r)/A) in 
each of the three regions. These expansions are shown below: 



ui 



u d 



3p 2 



2n 2 A 2 (p 2 +773-7/2 
3 



-D m (A)f l0 (p) 



1 + 



d 1J0 (A) Up) 



2tt 2 A 2 
3 

2tt 2 A 2 



D ldo (A)fdo(p) 
Diho(p, A)/«)(p) 



Didijp) fdijp) Didijp) fdijp) 
D ld0 (A) faoip) D ld0 (A) f m {p) 

hxip) 



1 + 



/ao(p) 



(5.119) 
,(5.120) 
(5.121) 



Note that the last term in contains the product of two first order corrections. Only 
other first order corrections resulting from this product should be kept. 

By using dimensional analysis, we can determine the pieces contributing to X and 
those contributing to X x . When integrated, the leading term in u\ is Oij)/ A) and thus 
becomes part of the equation for X\. The other terms in u\ may be dropped because 
they are 0(rf / A 2 ) or higher. These statements also apply to u d . Therefore, only the 
leading terms of U\ and u d need be considered. Together they provide a contribution 
of 



dp 



3p 2 



27T 2 A 2 (p 2 +T)l -r)l 



-D m {A)f m {p) 



2vr 2 A 2 



Dido(A) fdo(p) 



(5.122) 
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Because D ld0 (A) = Du (A), we can simplify the equation further to obtain 



3 f°° 

which we know from Eq. 1)5.81)1 is zero. 

The only contributions to X must then come from Uh- The first term is 0(1), 
making the second 0(i]/A). These two terms provide us with the final equations: 

3 r°° 

J o = 7T2K2 d P D lh0 (p,A)f h0 (p,A), (5.124) 

Z7T A JO 

3 r°° 

Jl = TTaI / d V D im { V ,A)f hl {r] 2 ^p,A). (5.125) 

With these equations, we now have everything needed to convert values of 5 to values 
ofG 3 . 



•pr 

p 2 +r] 2_ ^2 



fioip) ~ fdo(p) 



(5.123) 



5.5 Quick Reference 



All previously derived equations are provided here for convenience. Related functions 

are grouped together. 

Pseudo-Wavefunctions: 



fm(V2,V3,p) 



fdoip) 

fdlip) 



ATT 2 D l0 (r] 2 ,r]3,p) Jo 



dq 



q (p 2 + - rj 2 ) 

P (<? 2 + vi~ vi) 



( rj + 2p 2 + 2q 2 + 2pq\ 

ln iiWW^J /,ofeM) 

~q(p 2 + vi -vl) 



A7i 2 D m (r] 2 ,ri 3 ,p) Jo 
x In 



dq 



p{q 2 + rg- rg) 



7]l + 2p 2 + 2q 2 + 2pq 
r)i + 2p 2 + 2q 2 - 2pq 

Dii(rj2,rj3,p,A) 
D l0 (r}2,r) 3 ,p) 



fio(v2,m,p) 



p 



50 dq ^ I p 2 + q 2 + pq ^ 



ATT 2 D d0 (p) Jo q \p 2 + q 2 - pq 



P 



D dl (rj 2 ,p,A) 



(5.126) 



(5.127) 
(5.128) 



An 2 D d0 (p) \ D d0 (p) 

x(U(q) + U(q)) 



\ f°° dq 


In 


( p 2 + q 2 


+ pq\ 


I Jo q 




\^p2 _|_ g2 


-pq) 



(5.129) 
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/ho(p,A) 



P 



2p 2 + 2g 2 + 2pgz 



2n 2 D h0 (p, A) 

_ 3frj + 8fr 2 + 16 /%feA)\ fV \ 



fhl(V2,V3,P, A ) 



P 



2n 2 D m {p, A) 
xf hl (rj 2 ,r] 3 ,q,A) 



rOO rl 

^jJo dq L 



i C/ 2 (g + \p) U 2 (p + §<f) 

<7 : - 



S 

Analytic Solutions: 

fdo(p) 

/di(?72,P,A) 



p 2 D lh0 (p,A) 



D h0 (p, A) 
5i G 



2p2 + 2 g2 + 2 pqz 
Ao(p,A) 



k 2 G 2 #D m {p,K) 



2,1 



G 2 ,o_ 



A cos ( s In 



cos s In — + 9, 



P 



A 



p 



+ A w — cos j^s ln 
+ A dl0 ^cos (s hi 



+ A dll ^ln 



+ e hl 

+ #dl0 

cos [ s In ^ 



+ 0. 



D Functions: 

Dio(r]2,V3,P) 
Dn(r]2,V3,P, A) 
D d0 {p) 
D dl (rj 2 ,p,A) 

D h0 (p,A) 



| p 2 + ^2 _ ^ 



8^71 

fe| + 8h 2 - 16 ( Ip^rrj-r^ 
64^3/2 V A 



16tt 



V2 



+ 



3(h 2 2 + 8h 2 - 16) p 2 



(5.130) 



256tt 2 A 4 

+ 3^h 2 Ap ((8 + A 2 - 3h 2 p 2 ) 



(5.131) 



(5.132) 



(5.133) 



(5.134) 
(5.135) 
(5.136) 
(5.137) 



8^71 128^71-3/2 A 



(5.138) 
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D hl (r] 2 ) = 
Di Functions: 

D m (A) 
Dm(r] 3 ,p) 

D ld0 (A) 

D ldl (p) 

Di h0 (p,A) 



A (6 (hi + 4h 3 + 12) + h 2 (5h 2 + 12h 3 + 12)) 
576 v / 3tt 3 / 2 



(5.139) 



7T 



A (6 (hi + Ah 3 + 12) + h 2 (5hl + 12/13 + 12)) 



576^7^/2 



7T 



2 

I roc 



P 



47T 2 JO 

D lhl (p) = 
Two-Body Coupling: 

G 2 
G 2 fl 

G 2 ,i 

Three-Body Coupling: 

G 3 = A 4 g 3 
G 3) o z 



dq q 2 I dz 



U2(q+ 1 2 P)UMUMU^P + q) 
2p 2 + 2q 2 + 2pqz 



Ag 2 

128^2 7r 3 / 2 
3h 2 2 + 8h 2 + 16 

2048^2 tt 2 
(3h% + 8h 2 + 16) 2 VA 



5n 



G 2 , lo - D 2 fld /A 4 



G 3: i 



8i r 1 



D 



G 2fi 1 + G 2 ,i Xo — ^ #i ) G\ 



J Functions: 



dpD lh0 (p,A) fho(p,A) 

2tt 2 A 2 7o dp - Dlh °( p ' A )^ ll ( ?te ' 773 ' p ' A ) 
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2tt 2 A 2 7o 
3 f°° 



(5.140) 
(5.141) 

(5.142) 
(5.143) 

(5.144) 
(5.145) 

(5.146) 
(5.147) 

(5.148) 

(5.149) 
(5.150) 

(5.151) 

(5.152) 
(5.153) 



Functions: 



= 8^ Jo Pdp i qdq L dZ VWT2^i (5 - 154) 
A>,i = (5.155) 
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CHAPTER 6 



SOLUTION METHODOLOGY 



Once the integral equations have been derived, an approach for solving them must 
be decided upon. The approach used here is a combination of analytical and numerical 
methods. By analyzing the equations, we discover certain properties that help make 
the numerical calculations more efficient. This results in higher accuracy and lower 
computational costs. 

This chapter begins with a brief discussion of the steps needed to solve the integral 
equations. We point out quantities that one might compute and what information is 
necessary to carry out that computation. We also attempt to illustrate the relation- 
ships between the parameters in the integral equations. 

After this discussion, details of the numerical solution technique are given. The 
method involves discretizing the integral equations with an exponential spacing of 
points to significantly improve convergence. Ways of reducing the number of necessary 
points by decreasing the range of integration are also discussed. The final discretized 
equations are exhibited. 

Our analysis will be limited to the set of leading order integral equations. These 
are the equations relevant to the A — > oo limit. The equations for the first order 
corrections must be solved if one is interested in the cutoff-dependent behavior. Al- 
though we do not attempt to numerically solve for these corrections, the equations 
have been derived and all that remains is implementing a program to solve them. 
Many, if not all, of the issues discussed in this chapter are relevant to their solution. 

6.1 General Approach 

One of the leading order equations has already been solved, namely the one for fdo(p)- 
The solution has been determined to be of the form A cos(s m(p/A) + 9). Because 
any solution to the integral equations can be rescaled to give another equally valid 
solution, we choose to always scale our solutions so that A = 1. 
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The solutions for fio(p) and fhoip) are inherently linked by the solution for fdoip)- 
For p 3> (772, 773), the function f i0 must approach f d0 , and for p <C A, must approach 
fdo- In essence, these requirements act as boundary conditions for our solutions. 

Consider the case where the values of B2 and B3 are known for a given cutoff 
A. With just the energy values, we have enough information to solve the equation 
for fio(p). Once a solution is found, we can examine its behavior at large momen- 
tum. It must have a cosine behavior, and we can use this fact to fit the function 
cos(so ln(p/A) + 9) to the solution in this region. This allows the determination of 
the phase, which we will call 9i . 

For the function fhoip), we see from the integral equation that the solution depends 
upon two parameters, A and S. Since A is known, we can choose a value for 5 and 
then solve the equation. The result is a solution for fho that, for very small values of 
p, has a cosine behavior like fdo- Call the phase of this solution 6ho- 

Our boundary condition can then be stated more precisely as 



for some positive integer n. 1 This shows that the value chosen for 5 cannot be 
arbitrary as it must be a value that allows (j6.1J) to be satisfied. We see that choosing 
the bound-state energies determines a unique value of 9, which in turn determines a 
unique value for 5. The value of G3 can then be computed from 5. 

Now suppose that we would like to compute the other three-body bound-state 
energies in the spectrum. The three-body coupling has been fixed by our initial 
data, which is equivalent to saying that the phase has been fixed. The problem then 
becomes one of finding other values for B% that will yield a solution for fio with 
the same asymptotic phase 9iq. Thus we find that we can calculate a bound-state 
spectrum by solving the equation for fio and matching phases, never once needing to 
calculate fho- 

All other calculations are simply a matter of relating G3 and B3 via the interme- 
diate phase 9. 



Because a closed-form solution for either fio or f^o is unknown, we must resort to 
numerical calculations of these functions as well as any energies or couplings. A 
common method for solving an integral equation involves changing the integration 
into a sum over discrete points. The integral equation then becomes a matrix equation 
easily solved by standard methods. 

The reason why this is ir and not 2ir has to do with the equation for fig. For any low- momentum 
solution with a phase of we can multiply the function by -1 to obtain a solution with phase 8 + it. 
The phases are only unique up to a difference of tt. 



9 10 = 9 h0 ± mr, 




6.2 Numerical Methods 
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While this may appear straightforward, there are several practical issues to con- 
sider. For example, new limits on the integral equation must be determined. It is 
impossible to numerically integrate to infinity, so a suitable upper bound must be cho- 
sen. In our case, we will use a logarithmic integration range so a new lower bound to 
replace must also be determined. How the discrete points are chosen must be care- 
fully considered. Should we use the basic Simpson's Rule for numerical integration, 
or perhaps a more complicated quadrature method? 

Each choice is a compromise between accuracy and size. For example, we could 
choose an upper limit that is extremely large (ensuring that we are "close" to infinity) 
and discrete points that are closely spaced (minimizing errors in the sum). The trade- 
off is that using more points requires using a larger matrix. Using N discrete points 
will result in an N x N matrix with iV 2 elements. If all numbers are double precision 
decimals, even 8000 points would be enough to overwhelm a computer with 512 MB 
of memory. This does not even take into account the time needed to process such a 
matrix. Obviously, the goal is to obtain the desired accuracy with a minimal number 
of points. We will discuss a few methods that drastically reduce the number of points 
we need. 

In the following sections, we will assume that our goal is about 12 digits of accu- 
racy. This high accuracy may not be necessary for most leading order calculations, 
but it is essential when including the first order corrections. Besides directly ob- 
taining the equations for the A — > oo limit, one principle reason for expanding the 
three-body equation in powers of rj/A is to analyze the cutoff dependence. If we are 
attempting to study this behavior, we must be certain that our numerical errors are 
not larger than the corrections being studied, otherwise there is no way to distinguish 
the small corrections from the numerical "noise." By laying the groundwork for high 
accuracy in the leading order calculation, higher order calculations should be easier 
to implement. 



6.2.1 Transition to Mid-Momentum Function 

One way to limit the size of the matrix is to limit the range over which the function 
must be integrated. We know that fi and fho approach f d0 for p ^> rj and p C A 
respectively. This allows us to replace either function with cos (sq ln(p/A) + 9) in the 
appropriate range. The point at which we can make the switch is determined by the 
accuracy we desire. These limits are derived for the case of 12 digits of accuracy. 
For p^$> T), the equation for fi can be written as 



fio(V2,Vs,p) 



4(i + vi/p 2 -vi/p 2 



VSn (y/l + (277 3 2 )/(3p2) _ J ° 1 



X 



In 



' p 2 + q 2 + pq s 



•nl 



vi 



dq 



vi 



p 2 + q 2 



m, 



2p 2 + 2q 2 + 2pq 2p 2 + 2q 2 - 2pq 



fio(m,m,Q)- (6-2) 



62 



Here we have treated rj/p as a small quantity and perturbatively expanded all factors. 
As long as both i] 2 /p and 773/p are less than 10~ 12 , this equation will match the one 
for fdo to 12 digits. Of course, 773 must be greater than r] 2 since we are considering 
only stable bound states. This means that p ~ 10 12 1] 3 sets the limit above which fi 
can be replaced by f d0 . In practice, we must have enough data points above this limit 
to ensure that our cosine fit is accurate to 12 digits also. Therefore, we will use an 
actual limit of p = 10 15 773. 

Similarly, we can expand the equation for fh in the region p < A. Notice that 
in this region D lhQ (p, A) approaches Di d0 (A) which is proportional to A, and the 
function D h0 (p, A) becomes equal to D d0 (p) = \^3p/l6n. In the f h0 integral equation, 
the momentum-dependent part of the three-body interaction becomes 

(6.3) 

Since the leading order mid-momentum equation has no three-body interaction term, 
the above term will equal zero to 12 digits if we choose p ~ 1CT 12 A. 
The integral part for f h0 looks like 

' dq / dz \ 9 2 V ? ,o 2 J ho(q,A). (6.4) 




2tt 2 D hQ (p, A) Jo J-i 2p 2 + 2q 2 + 2pqz 

We have already stated that D h0 (p,A) approaches D d0 (p), but more importantly, it 
approaches like 



D h0 (p,A) p -^D d0 (p) 



x f V3(h 2 2 + 8h 2 -16) \ p 
\ S v '2^ ~ ) A 



(6.5) 



It will therefore equal D d0 to 12 digits if p ~ 1CT 12 A. This limit also applies to the 
integrand itself, so we may replace fho with f d0 for values of p less than (9(10 -12 A). 
In practice however, we use a limit of p = 1CT 17 A to ensure that we have enough 
points below this region to fit the cosine behavior. 



6.2.2 New Integration Limits 

For the case of fho, we have limited the range of integration to be 10~ 17 A to 00. 
(Below this range, we use f d o-) Naturally, we cannot integrate to infinity and instead 
must find a new limit to replace it. Let us call this limit A. We choose A such that 

v 2 r°° r 1 U 2 (q + \p] U 2 (p + \q] 

2^L(M) L d " L d * y + ^ A) < 0(10 >' (6 ' 6) 

The exponentials in U 2 suggest that the integrand should die off quickly, allowing 
us to make an initial guess for A using 

g— A 2 /A 2 = 1Q -12 (6J) 
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This gives an initial value of A ~ 5.25 A. However, the double integral makes the 
analysis harder since we cannot determine the exact behavior. We must resort to 
numerical computation, and some sample calculations reveal that a limit of A = 10 A 
is sufficient for our purposes. 

From 10 _17 A down to 0, fho is replaced by fao- We would like to replace the 
limit with a larger value that still maintains our desired accuracy. Even though we 
know the analytic solution for /^ , narrowing the range of integration will reduce our 
computational effort. Call this new lower limit e, which is chosen so that 



2ir 2 D h0 (p 



i U 2 (q + \p) U 2 (p 



2p 2 + 2q 2 + 2pqz 



f d0 (q,A)<O(W 



-12\ 



(6.8) 



We assume that e 10 17 A < p. To within 12 digits of accuracy, the z integration 
can be written as 



i U2 (q + \p) U2 (p+ \ q 
-1 Z 2p 2 + 2q 2 + 2pqz 



1 dz U 2 {p/2)U 2 {p) 
1 2p 2 + 2pqz 



1 



= -C/ 2 (p/2)f/ 2 (p)ln 



p + pq 



m, 



(6.9) 



Since q p, the logarithm can be approximated as 2q/p. Our constraint for e now 
becomes 

2^( P ,A) (P72 > lh ® jf dqUq) < ° (10 " 2) ' <6 ' 10) 

The values of p are of the same order as A, so we expect the value of U 2 to be 0(1). 
The function f^o is also 0(1), so it is simply replaced by 1 in this approximation. This 
leaves an integral with a value of e. When combined with D h0 (p, A) ~ 0(p), we find 
that e/p < 0(1O~ 12 ). The smallest value for p is 10~ 17 A, implying that e ~ 10^ 29 A. 
In practice, this value is sufficient for 12 digits of accuracy. 

Having replaced the limits for f^o, we move on to fio- Again, we must find a finite 
upper limit to substitute for infinity. For p > 10 15 773, fio is replaced by fdo- This new 
limit, A, is determined by the condition 
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fdo(q) 
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(6.11) 
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where we have used the fact that q ^> p ~ 773. For values of p much larger than 
773, the term in Eq. (J6.11j) is proportional to p/X. The largest value p can obtain is 
10 15 773, implying A = 10 27 r/ 3 . Numerical calculations verify that this limit is sufficient 
to assure the desired accuracy. 

Finally, we must replace the lower limit for fio with a non-zero value e that we 
assume to be much smaller than 7/3. Our requirement is that 

(p 2 + rj - 77I) dqq\n ( V * + 2 ^ + 2PQ 

4n 2 p(r] 2 - r)%)D l0 (r) 2 , 773, p) Jo \T)i + 2p 2 - 2pq 

(6.12) 

Expanding the logarithm yields 

2 w tl\\Tn ( T f e dqq 2 HV2,V3,q) < 0(1O- 12 ). (6.13) 

v 2 m - m)\m + 2p 2 )Ao(^2,^ 3 ,p) 

For small values of q, we will find that fio is approximately constant and of 0(1). 
Therefore, the integral is roughly equal to e 3 . If p <C 773 or p ~ 773, then Eq. (j6.13|) 
is 0(e 3 /r] 3 ). If p 3> 773, then it is 0(e 3 /pr] 3 ) < 0(e 3 /rf). This seems to imply that a 
value of e ~ 10 — 4 773 is adequate. 

Unfortunately, using this value of e will result in poor accuracy when 772 — 
This is a result of the (rj 2 — r/|) term in the denominator of Eq. ()6.13|) . Originally, 
this was our approximation to the term (q 2 + rj 2 — 77!) in the integral equation. When 
the energies are nearly equal, our approximation needs to be q 2 . The condition on e 
should now become 

/ 2 2 2 \ 

V 2 I 9 + 2 ? 3 n7 V f e dqf lQ (V2,m,q) < 0{lQ- 12 ). (6.14) 

ir 2 (r]i + 2p 2 )D l0 (T] 2 , 773,79) Jo 

The integral is roughly equal to e, and the entire term is 0(e/r]). This means we must 
use the lower limit e = 10~ 12 773. 

These derivations are general enough to determine the appropriate limits for other 
cases of desired accuracy. Keep in mind however that the limits must always be tested 
numerically to ensure that they are indeed sufficient. 

6.2.3 Discrete Point Spacing 

Now that we have stricter limits in place, we must choose the discrete values of 
q within these limits at which to evaluate our functions. We employ discrete points 
that are equally spaced on a logarithmic scale. There are two main reasons for making 
this choice. 

First, we have already seen that f d0 = cos (s ln(p/A) +9) is periodic on a loga- 
rithmic scale. In fact, other functions have similar behavior, including G3. It makes 
sense that a logarithmic spacing is better suited to capturing the behavior of the 
system. 



fi (V2,V3,q)<O(W 



-12\ 
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Second, by choosing points in this manner, an integration of q from to oo becomes 
an integration of ln(g) from — oo to oo. By equally spacing points on this log scale, we 
can achieve convergence that improves exponentially 2 with the spacing, as opposed 
to the more typical power law convergence. This drastically reduces the number of 
points needed to achieve our desired accuracy. For instance, to cover the range of fio 
from 10~ 12 773 to 10 15 773, and maintain 12 digits of accuracy, we can space our points 
by Pn+i = Vn,e°' 2 - This requires only about 300 points. 



6.3 Discretized Equations 

The discretized integral equations are shown below. Here, A represents the momen- 
tum spacing, and the momentum values are related by the equation p n+ \ = p n e A . 
The identity matrix is represented by I nm , and the ceiling function represented by 
\x] returns an integer value n such that (n — 1) < x < n. 

Low-Momentum: The value of n for p n ranges from to N m id, while the value 
of m for q m ranges from to N max . These are defined as: p = 10~ 12 ^ 3 , N, 



mid 



[ln(10 27 )/A], and N max = [ln(10 39 )/A] . 

{M nm - I nm ) fioipm) = b n (6.15) 



m=0 



^ 2 D m (r] 2 ,7]^p n )p n (q 2 m + ril-ril) \ + 2p 2 n + 2q 2 m + 2p n q m J 



^{pl + vl-vl 



An 2 p n D l0 (r] 2 , rj 3 , p n ) m= ^. d+1 



ln 



m + 2 Pn + 2 <lm + 2 P»gr, 

Vi + 2 Pl + 2q 2 a + 2p n q ri 
xcos Lin (jfj+ej 



(6.17) 



High-Momentum: The value of n for p n ranges from to N max , while the 
value of m for q m ranges from —N min to N max . These are defined as: p = 10~ 17 A, 
N max = [ln(10 18 )/A], and N mm = [ln(10 12 )/A] . 

^ (M nm - I nm ) fmiPrn) = a n + 5 b n (6.18) 

m=0 

2-K z JJ h0 {p n ,A) 

P ' A £ q m K( Pn , q m ) cos (s In ( + o) (6.20) 



2n 2 D h0 (Pm A) m ^ N _ *™ ^ ™' V- V A 



2 For a basic mathematical proof of this claim, see Appendix [B] Appendix O has a short error 
analysis that numerically shows the exponential convergence. 
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CHAPTER 7 



ANALYTIC AND NUMERICAL RESULTS 



A computer program to solve the discretized integral equations can be imple- 
mented with the help of some basic numerical algorithms. Once completed, it can be 
used to generate highly accurate results and analyze the three-body system. Nonethe- 
less, attempts to extract analytic results should not be overlooked. Even some of the 
most general properties of the integral equations allow us to draw conclusions about 
the behavior of the system. 

We begin this chapter with analytic results obtained from studying the leading 
order integral equations. These results include statements about the cutoff depen- 
dence of bound-state energies and the phase for fdo- A proof for the cyclic behavior 
of 5q is given, which is then used to infer similar behavior for G3. 

Following the analytic results is a section containing numerical solutions to the 
integral equations. Here we examine behavior that cannot be determined analytically. 
Solutions for the functions fi and fho are shown, and the cutoff dependence of G 3 is 
calculated. Some relations that are proven analytically are also verified numerically. 

7.1 Analytic Results 

In Chapter we saw that Eq. f)5. 126j) contains no A dependence, making the function 
fio(p) cutoff-independent. However, for values of p ^> 773, the solution must behave 
like fdo(p) = cos (s In (p/A) + 9), which explicitly has the cutoff in it. The only way 
these two equations can be reconciled is if 9 contains cutoff dependence of the form 
So ln(A). Any remaining part of 9 must be a function of r] 2 and 773. Therefore, we 
shall write 

9 = s In (A/773) +HV2/V3), (7.1) 

where 9 is a dimensionless function of the ratio 772/773. This relation holds for any 
values of 772 and 773, including all 773 values corresponding to multiple bound states with 
the same 772. In the next chapter, we shall see that 9 is closely related to Efimov's 
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universal function [T^]. Using 773 in the ratio with A is a matter of choice. Any 
quantity composed of r] 2 and 773 with the same dimension as A would work just as 
well; however, we need to allow r] 2 = 0, so r] 2 alone is a poor choice. 

While Eq. (17. 1(1 is quite simple, it has many interesting consequences. As we 
mentioned earlier, two different three-body bound states in the same spectrum must 
have the same phase 9. Suppose that we choose some fixed values for 772 and 773 and 
make them cutoff-independent by choosing the appropriate A dependence for G 2 and 
G3. This might be desirable if we are trying to match those energies to experimental 
data. The phase for the solution in this case is 

6 A = s In (A/773) +HV2/V3), (7.2) 

for some given cutoff A. If the cutoff is changed, G 2 , G3, and 6\ will all change with 
it, but 772, 773, and #(772/773) will not. 

Imagine now that we find a second three-body bound-state solution with the same 
phase. Let us call its energy 773. The phase for this solution is 

0A = soMA/7fe) + 0(W*fe), (7-3) 
which must be equal to #a by assumption. This results in the relation 

s In (7/3/ 773) = (772/773) - 6 (772/773) • (7.4) 
If the cutoff is now changed to a new value A', the original data gives a phase of 

0A' = soMA'/»fe) + 0(V%). (7-5) 
The question is whether 773 is still a valid solution. Its new phase is 

K > = s In (A'/t/g) + 6 (773/773) 

= s In (A'/ 773) + [s In (773/ 773) + 9 (r/ 2 / 773) 

= s In (A'/r/s) + ^ (7/2/773) 

= A ,. (7.6) 

Since the phases still match, 773 is still a bound-state solution. This remains true for 
any cutoff, implying that 773 is also cutoff-independent like 773. Of course, the same 
statement applies to any other three-body bound state making the entire spectrum 
completely independent of A. Such behavior should come as no surprise since the 
leading order equations represent the A — > 00 limit. Keep in mind that this is true 
for any other physical quantity, but does not apply to the couplings. Obviously, it 
is the cutoff dependence of the couplings that enables the bound states to be cutoff 
independent. We have shown that a single three-body contact interaction allows us 
to renormalize the entire three-body bound-state spectrum. 

Just as fio has no dependence on A, neither does it have any dependence on h 2 . 
As a consequence, h 2 has no effect at leading order on the binding energies or other 
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physical quantities. It does appear in the equation for fn however, showing that it 
is needed when considering first order corrections. Because we are working only to 
leading order, we shall use h 2 = in our calculations from now on unless otherwise 
specified. 

We now turn to the equation for fhoip)' 

V 2 r 1 u 2 (q + hp) U 2 (p + lq) 

/ftofoA = P / dq dz \ 2 " J 2 M \ y Ao( 9 ,A) 

27i z D h0 (p, A) Jo J-i 2p l + 2q l + 2pqz 

__ 3h 2 + 8h 2 + 16 ( D lh0 (p,A) \ 

° 128^^/2 {D h0 (p,A)){A*)- { '-' } 

The dependence on p and A has been explicitly shown. Since fto is dimensionless, 
only the ratio p/A can occur in the function. Furthermore, there is a dependence 
upon h 2 and 9. The h 2 dependence comes from its appearance in the function U 2 (p), 
either directly in the integral or indirectly in Dho and -Di^o- The 9 dependence is a 
result of fho approaching fdo for p ^ A. 

Assuming that A is held fixed, choosing a value for 5q will uniquely determine the 
phase. Thus, we can view the phase as a function of the coupling, 9(S ). Conversely, 
choosing a phase determines the coupling, so it is equally valid to treat the coupling 
as a function of the phase, S (9). The coupling 5 can also have a dependence upon h 2 
but not upon A. The reason is that 5o is dimensionless and there is no other quantity 
available to form a dimensionless ratio with A. 

With this in mind, consider a solution to Eq. ()7.7j) with a coupling of £o(#) and a 
function fho(p) that behaves asymptotically as cos (so In (p/A) + 9) for small p. For 
a phase of 9 + 2tt, the cosine behavior will remain unchanged. We may therefore 
conclude that 

5 (9 + 2n) = 5 (9), (7.8) 

showing that the coupling exhibits periodic behavior. If we apply the operator 
—d 2 /d9 2 to both sides of (j7.7J) . we find once again that the asymptotic cosine be- 
havior is the same. This implies 

- H = M«). (7.9) 

The conclusion to be drawn is that the coupling may be written as a cosine function 
with some amplitude A and phase <fi. These two parameters will contain any h 2 
dependence that 9 may possess, so we shall write 

5 (8) =A{h 2 ) cos(# + 0(/i 2 )). (7.10) 

The equation relating <5 and G3, 

G3 ' 0= G 2 , J -\oVA 4 ' (7 - 11} 
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shows that whenever Sq is zero, so is G3. 1 Since two adjacent zeros of Sq occur when 
the cutoffs satisfy 



the adjacent zeros of G 3 (A) should be spaced by a cutoff factor of A'/A = exp(7r/s ) — 
22.69438259536. This suggests that G 3 (A) may also possess cyclic behavior, but this 
must ultimately be verified numerically. 

7.2 Numerical Results 

We begin our numerical investigation by considering solutions for the functions fa and 
fho which lead to an approximation for the complete function f(p). Next, the cutoff 
independence of the three-body spectrum is verified, followed by an analysis of the 
coupling constants. These numerical results are confirmed by Wilson's calculations 

B5|. 

7.2.1 Low- and High-Momentum Functions 

Figure ITTT1 shows the numerical solution for fm(p) in the case of B 2 = 0.1, B3 = 1.0, 
and A = 10 8 . Notice that it is constant for small values of momentum and then takes 
on the cosine behavior as p becomes large. This behavior is typical of all solutions 
for fi . For reference, we have included numerical solutions for the additional cases 
of B 2 = 0.0,0.5, and 1.0 in Figs. 17.21 17.31 and 17.41 respectively Since the behavior of 
the solution to the low-momentum equation is determined by the ratio B2/B3, these 
figures provide a representative sample of all solutions. We can also study the effects 
of changing the cutoff while keeping the energies fixed. Equation (|7.1|) shows that 
changing the value of A will only result in a shift in the value of 9, which is exactly 
the behavior seen in Fig. 17.51 Note that we follow the convention of restricting all 
phases to be within the range to 2n, which explains why one of the phase values is 
close to zero. 

Several bound states may exist for the same values of B 2 and A. In the case of 
B 2 = 1.0 and B 3 = 1.0 for A = 10 8 , there exist states of energy B 3 = 6.7502901502599 
and B3 = 1406.130393204. The solutions of fio for these energies are shown in Fig. 17.61 
All of the functions are very similar. The only real difference is in the length of the 
initial "plateau." As the bound-state energies become larger, the flat region becomes 
longer, joining the cosine part near later peaks. An additional illustration of this 
behavior is shown in Fig. 17.71 using different bound-state energies. 

Using the parameters 9 = 5.684386276089572 and A = 10 8 , we find the numerical 
solution for the function fho{p) shown in Fig. 17.81 This is the high-momentum function 

lr Tlie only problem that might arise is if 2q — in such a way that the ratio is non-zero, but we 
have found no parameters for which this is true. 




(7.12) 
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corresponding to the low-momentum function in Fig. 17.11 Figures 17.91 to I7TTT1 display 
the high-momentum counterparts of Figs. 17.21 to 17.41 For all fho solutions, we see the 
cosine behavior for p <C A and a suppression of large momentum values when p > A. 
This suppression is an effect of the gaussian behavior of U 2 . The behavior of fho is 
determined by the phase. Different values of the cutoff for the same phase should 
result in the same function. An example of this is shown in Fig. 17.121 where 9 = 3.0 
for A = 10 5 , 10 8 , and 10 11 . Notice that all three solutions coincide. 

Combining f [0 , f d0 and fho, we can see the overall leading order behavior for the 
function f(p). Plots of // + fho — fdo are shown in Figs. 17. 131 to 17. 161 for the four low- 
and high-momentum solutions previously considered. These solutions span several 
values of B 2 /B 3 and show the typical behavior of any leading order solution. Figure 
17. 171 shows the three individual components for B 2 = 0.1 on the same plot. From this, 
it easy to see how cancellations occur in each region making it clear why fi + fh — fd 
is a good approximation to / everywhere. 

7.2.2 Three-Body Binding Energies 

In Sec. I7.1l we proved that the three-body bound-state spectrum is cutoff independent 
to leading order. This is verified numerically in Fig. 17.181 Here we have chosen B 2 = 
and let G3 change with A so that the state B 3 = 1.0 is held constant. Two other 
states, one shallower and one deeper, are calculated as the cutoff changes. Since very 
small fluctuations are impossible to see in the plot, we have included Table l?~T| which 
shows the calculated energies for several values of A. Using Efimov's result that the 
ratio of adjacent binding energies is exp(27r/s ) when B 2 = [EJ|2I|, the relative 
error for each calculation can be determined and is also given in the table. This 
illustrates that each energy is cutoff-independent to about 12 digits and also matches 
the true value to the same accuracy. As an additional example, Table 17.21 shows the 
case of B 2 = B 3 = 1.0 and considers the next two deeper states as A changes. The 
binding energies are approximately 6.75029 and 1406.13. These energies have been 
previously calculated by Braaten, Hammer, and Kusunoki [23] using a method that 
gives at most two digits of numerical precision. Their results are 6.8 and 1.4 x 10 3 , 
which match the results given here to within relative errors of O(10~ 3 ). Note that 
our values for these states are independent of A to about 12 digits. 

We have shown that a three-body interaction is capable of fixing a three-body 
binding energy. However, this does not prove that such an interaction is necessary. 
To investigate this claim, let us assume that G3 = for any value of A, so that the 
three-body interaction is removed from our equations. Using B 2 = 0, we calculate 
two values of -B3 for a cutoff of A = 10 5 . The energy of each state is tracked as the 
value of A changes. The results are shown in Fig. 17.191 
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A 


B 3 (Shallow) 


Error 


B 3 (Deep) 


Error 


100000.00000000 
738905.60989306 
5459815.0033144 
109663315.84284 
3631550267.4246 


0.0019416156131338 
0.0019416156131358 
0.0019416156131435 
0.0019416156131358 
0.0019416156131435 


6.7e-13 
3.2e-13 
4.3e-12 
3.2e-13 
4.3e-12 


515.03500138461 
515.03500138403 
515.03500138287 
515.03500138520 
515.03500138520 


5.2e-13 
1.6e-12 
3.9e-12 
6.0e-13 
6.0e-13 



Table 7.1: Binding energies of the next shallowest and next deepest 3-body bound 
states for B2 = 0.0, B% = 1.0 and various cutoffs. Several data points from Fig. 17.181 
are shown along with the relative error of each result. 



A 


5 3 #1 


B3 #2 


100000.00000000 
738905.60989306 
2008553.6923187 
14841315.910257 
298095798.70417 
5987414171.5197 
44241339200.892 


6.750290150257678 
6.750290150257678 
6.750290150255419 
6.750290150268966 
6.750290150257678 
6.750290150259935 
6.750290150257678 


1406.13039320296 
1406.13039320593 
1406.13039320345 
1406.13039320345 
1406.13039320593 
1406.13039320345 
1406.13039320296 



Table 7.2: Binding energies of the two next deeper 3-body bound states for B2 = 
-B3 = 1.0 and various cutoffs. 
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Both binding energies diverge like A 2 . This is easily seen using equation (|7.1|) . 
Since G3 = 0, the phase must be fixed. This means 



from which we can easily derive the relation B 3 oc A 2 . This divergent behavior is an 
example of the unbounded three-body spectrum predicted by L. H. Thomas for short- 
ranged interactions [21 J. While this alone is not a proof that a three-body contact 
interaction is required, it should be possible to show that any other local interaction 
does not enter the leading order equation. This would imply G 3 is the only coupling 
capable of making any binding energy independent of the cutoff. 

7.2.3 Couplings 

Equation ()7.10|) shows that the coupling S should have a cosine dependence on the 
phase 9, which is defined by fao = cos(s ln(p/A) + 9). Numerical data for 5 as a 
function of the phase is shown in Fig. I7.2U1 along with a best-fit cosine curve to verify 
the behavior. This behavior is independent of B2 and B3. Figure 17.211 proves this by 
plotting the same relation for different energies. Keep in mind that /i2 can affect the 
amplitude and/or phase shift of the cosine behavior. 

We suggested in Sec. 17.11 that this periodic behavior should carry over to G 3 . 
Figure 17.221 displays a plot of G 3 as a function of A. This data exhibits the limit- 
cycle behavior of the three-body coupling. As the cutoff increases, G 3 becomes larger 
and larger, eventually diverging to infinity. It then jumps to negative infinity and 
continues increasing. The data for 5q and G3 are combined in Fig. 17.231 to show that 
5q = implies G3 = 0, a claim also made in Sec. 17.11 Figure l7~24l uses two different 
sets of energies to show that the limit-cycle behavior is not dependent upon any 
specific bound-state values, but the positioning of the cycle is dependent upon the 
energies. 

If the scattering length a is large compared to the effective range, the zero angular 
momentum sector of the three-body problem with short-ranged interactions reduces 
to a one- dimensional problem with a potential of 1/r 2 in the region r e <ti r ^ a [26j. 
The 1/r 2 potential is studied in [27] where a short-range interaction is used to regulate 
the divergences from the long-range 1/r 2 interaction. While one might expect to find 
cutoff-dependent behavior similar to what we see, they instead find a monotonically 
increasing coupling with no divergences. They suggest that the limit-cycle behavior 
is not a universal aspect of the renormalization group flow but rather an artifact of 
the regularization method used. 

However, this cyclic behavior has been previously observed ^2] using a sharp 
cutoff that simply discards all momenta above A. Our method does not discard such 
momenta. It merely suppresses them by a gaussian factor. Yet the limit cycle still 
remains. This evidence seems to indicate that the limit-cycle behavior is universal, 
and the issue deserves further investigation. 




(7.13) 
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Another coupling constant we can study is hi. It is shown above that h 2 has no 
effect on the binding energies to leading order, but it may affect other quantities like 
G 3 and 5 . In fact, Eq. (|7.1U|) explicitly exhibits such dependence. In Figs. 17.251 and 
17.261 we plot the effect of changing hi upon So and G3 respectively. The range of h 2 
values is determined from Eq. ([2.34)1 . reproduced here for convenience: 

1 Z + (**U + (* + L*) =0 ( , 14) 



For very large cutoffs, the effective range as computed from hi takes the form 

J_ (hl + 8hi- 16). (7.15) 



2V27rA 

For illustrative purposes, we restrict ourselves to the case r e > so that the values of 
h 2 are limited to the range —4(1 + v2) < h 2 < —4(1 — y/2). This example is chosen 
simply to provide a window of hi values to concentrate on. 

The cosine behavior of Sq(6) should remain unchanged if h 2 is changed, but the 
amplitude and phase of the curve may differ. This is exactly what we find in Fig. 17.271 
for the cases hi = 0.0, —2.0, and —5.0. It appears as though all curves intersect at 
the same two points. However, we do not know if this can be explicitly proven from 
the leading order equations or what its significance may be. 

The general limit-cycle behavior of G3 also stays the same for hi = 0.0, —2.0, and 
—5.0. These curves are shown together in Fig. 17.281 Changing hi results in a shift 
and/or flattening of the limit cycles, but the periodic behavior is still apparent. 
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Figure 7.1: Numerical solution of for the case of £?2 = 0.1, B 3 = 1.0, and 

A = 10 8 . The dashed line is the best-fit cosine curve that matches the high-momentum 
behavior. 
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Figure 7.2: Numerical solution of fio(p) for the case of B 2 = 0.0, B 3 = 1.0, and 
A = 10 8 . The dashed line is the best-fit cosine curve that matches the high-momentum 
behavior. 
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Figure 7.3: Numerical solution of fm(p) for the case of B 2 = 0.5, B 3 = 1.0, and 
A = 10 8 . The dashed line is the best-fit cosine curve that matches the high-momentum 
behavior. 
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Figure 7.4: Numerical solution of fm(p) for the case of B 2 = 1.0, _B 3 = 1.0, and 
A = 10 8 . The dashed line is the best-fit cosine curve that matches the high-momentum 
behavior. 
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Figure 7.5: Numerical solution of fio(p) for the cutoffs A = 10 5 , 10 8 , and 10 11 with 
energies B2 = 0.1 and B3 = 1.0. The phases corresponding to these cutoffs are 
5.016726935003937, 5.684386276089572, and 0.06886030999563049 respectively. 
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Figure 7.6: Numerical solutions of fio for several bound states using B2 = 1.0 
and A = 10 8 . The three-body bound-state energies are 1.0, 6.7502901502599, and 
1406.1303932044. 
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Figure 7.7: Numerical solutions of fio for several bound states using B2 = 0.01 
and A = 10 5 . The three-body bound-state energies are 0.2, 63.033762419242, and 
31645.59444559555. 
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Figure 7.8: Numerical solution of f h0 (p) for the case 9 = 5.684386276089572 and 
A = 10 s . The dashed line is the best-fit cosine curve that matches the low-momentum 
behavior. 
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Figure 7.9: Numerical solution of f h0 (p) for the case 9 = 5.299886863265019 and 
A = 10 s . The dashed line is the best-fit cosine curve that matches the low-momentum 
behavior. 
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Figure 7.10: Numerical solution of f h0 (p) for the case 9 = 0.1031320389575355 and 
A = 10 s . The dashed line is the best-fit cosine curve that matches the low-momentum 
behavior. 
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Figure 7.11: Numerical solution of f h0 (p) for the case 6 = 1.681619491361092 and 
A = 10 8 . The dashed line is the best-fit cosine curve that matches the low-momentum 
behavior. 
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Figure 7.12: Numerical solution of fm{p) with 6 = 3.0 for cutoffs A = 10 5 , 10 8 , and 
10 11 . 
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Figure 7.13: Numerical solution of fo(p) for the case of B 2 = 0.1, -B 3 = 1.0, and 
A = 10 8 . 
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Figure 7.14: Numerical solution of fo(p) for the case of B 2 = 0.0, B 3 = 1.0, and 
A = 10 8 . 
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Figure 7.15: Numerical solution of f (p) for the case of B 2 = 0.5, B 3 = 1.0, and 
A = 10 8 . 
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Figure 7.16: Numerical solution of fo(p) for the case of B 2 = 1.0, B 3 = 1.0, and 
A = 10 8 . 
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Figure 7.17: Numerical solutions of fa, fdo, and fho for the case of B 2 = 0.1, B 3 = 
and A = 10 8 . 



92 



10000 



1000 



■xxxxxxxxxxxxxxxxxxx> 



100 



10 



DD 



0.1 



0.01 



0.001 



0.0001 




100000 



1e+09 



Figure 7.18: Binding energies of the next shallowest and next deepest 3-body bound 
states as a function of the cutoff A for B 2 = and B 3 = 1.0. 
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Figure 7.19: Binding energies of two three-body bound states for the case of B2 = 
and G3 = 0. Notice that without a three-body interaction, all the binding energies 
diverge as the cutoff is increased. 
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Figure 7.20: The coupling 5q as a function of the phase for the case B2 = 0.5 and 
B 3 = 1.0. The dashed line is the best-fit cosine curve. The data matches to within 
an error of about 10~ 10 . 
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Figure 7.22: The dimensionless three-body coupling G3 as a function of the cutoff A 
for B2 = 0.5 and B% = 1.0. The limit-cycle behavior of G3 is evident from the fact 
that it is periodic in ln(A). 
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Figure 7.23: Dependence of G3 and So on the cutoff A for the case of B 2 = 0.5 and 
i? 3 = 1.0. The values for 5 have been multiplied by a constant factor of 150 to make 
the cosine behavior visible. Both couplings are shown to emphasize that they equal 
zero at the same points. 
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Figure 7.24: G 3 as a function of the cutoff A for two different sets of B 2 and B 3 . Both 
curves have the same limit-cycle behavior, but the difference in energies causes one 
to be shifted relative to the other. 



99 



20 



15 



10 - 




-10 - 



-15 - 



_20 I 1 1 ± 1 1 1 

-10 -8 -6 -4-2 2 

h 2 



Figure 7.25: The ratio S (h 2 )/S (0) as a function of h 2 for 6 = 1.0 and A = 10 5 . 
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Figure 7.26: The ratio G 3 (h 2 )/G 3 (0) as a function of h 2 for 6 = 1.0 and A = 10 5 . 
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Figure 7.27: The coupling 5 as a function of 9 for several values of h 2 - The curves 
are all cosines with different amplitudes and phases. 
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Figure 7.28: G3 as a function of A for several values of hi- A change in hi results in 
a shift and/or flattening of the limit-cycle curve, but the general periodic form is the 
same. 
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CHAPTER 8 



EFIMOV'S FUNCTION 



In the 1970's, Vitaly Efimov investigated the behavior of three-body systems whose 
particles interacted via pair-wise "resonant" forces. A resonant force is one that 
supports a shallow two-body bound state or a virtual state close to the two-body 
threshold, resulting in a scattering length whose absolute value is much larger than 
the range of the interaction. This disparity in length scales allowed him to discover 
properties that were independent of the exact form of the two-body interaction. One 
of the most useful results is the relation of the entire three-body spectrum to a single 
universal function. 

This chapter begins with a brief introduction to Efimov's work. We outline the 
arguments used to derive Efimov's universal function and exhibit the discrete scaling 
symmetry it possesses. Some uses of his work in the fields of nuclear and atomic 
physics are discussed. Since our three-body problem satisfies the constraints used 
by Efimov, we use the leading order integral equations to numerically compute Efi- 
mov's function to high accuracy. We also analytically derive the relationship between 
Efimov's function and the function 9 defined in Eq. (J7.1j) . 

8 . 1 Background 

Efimov's work hinges upon the use of a resonant two-body force. Assuming a positive 
scattering length, this simply means that the potential associated with this force has 
a shallow two-body bound state. Since B 2 = h/ma 2 to leading order, this is the 
same as saying the potential has a large scattering length. Such a situation may arise 
either naturally (e.g., the deuteron) or artificially (e.g., alkali atoms by fine-tuning an 
external magnetic field). In any case, it will be assumed that the scattering length is 
much larger than the range of the interaction, which we call r . Efimov also assumed 
that there are no deeply-bound states in the two-body spectrum. 
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We begin with the adiabatic hyperspherical representation of the three-body equa- 
tion j26[ |2Ej • The hyperspherical radius R is defined by 

R 2 = \ (^2 + ^3 + 40, 

where = [ — r^- 1 . In the region ro -R a, the three-body Schrodinger equation 
for a zero angular momentum state takes the form 

h 2 \d 2 s 2 + 1/4" 



2m 



dR 2 + R 2 



f(R)=Ef(R). (8.2) 



Here, f(R) is the radial wavefunction and sq ~ 1.00624 is the same constant seen 
in previous chapters. The goal is to construct solutions for the bound state energies 
E = -B 3 . 

Since we are interested only in low-energy states (\E\ ~ I /ma 2 ), the requirement 
R ^ a implies that the 1/R 2 potential is much larger than E. Therefore, we can 
ignore the right-hand side of Eq. ()8.2|) . The resulting equation has a general solution 
of the form 

f(R) = VHR [A e is ° HHR) + B e~ [so ln( ^ } ] , (8.3) 

where H has dimensions of momentum, and A and B are dimensionless. For bound 
state solutions, these quantities must be functions of the bound-state energy B 3 and 
the two-body scattering length a. Let us define 



' mB% 1 



£ = — a.Tcta.n(aymB 3 /h). (8.5) 

This definition of £ is valid only for positive scattering lengths. For a < 0, we must use 
the definition £ = — a,Tcta,n(a\/mB 3 /h) — it. In either case though, the dimensionless 
quantities A and B can depend upon £. 

The radial wavefunction in the region ro C C a is a sum of incoming and 
outgoing hyperspherical waves. At short distances (R ~ Tq) and long distances (R ~ 
a) the wavefunction is more complicated, and it is the boundary conditions at these 
points that determine the actual bound-state energies. 

We first consider the constraints from unitarity in the short-distance region R ~ 
r . Because there are no deeply-bound two-body states, all the probability associated 
with the incoming hyperspherical wave must be reflected into the outgoing wave. This 
implies that A and B can only differ by a phase: B = — Ae ld . The wavefunction in 
the region ro <C R a is now 



f{R) oc VHR sin (so hi(HR) — 6/2) , (8-6) 
so we may express the phase as 

,f'(R) i y 



-Sq \n(HR) + arccot 



R 



f(R) 



•7) 
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The phase is determined by specifying the derivative R f > (Ro)/f(R ) at any point 
Ro satisfying r -C fl < a. We will simplify the form by writing it as 

= 2s o ln(#/A*). (8.8) 

The quantity A* is a complicated function of _R an d Rof'(Ro)/ f(Ro), the details of 
which are unimportant. What is important is that the derivative at Rq is ultimately 
determined by the short-distance behavior. Essentially, A* parameterizes the effects 
of all short-distance interactions without knowledge of the detailed forms. 

We now consider the constraints from unitarity in the long-distance region R ~ a. 
Because we are considering three-body bound states, all of the probability associated 
with the outgoing hyperspherical wave must be reflected into the incoming wave. This 
implies that A and B differ only by a phase: B = — Ae lA . The phase A is determined 
by the long-range parameters B 3 and a. We will write this phase as A(£) to explicitly 
show its dependence upon B 3 and a through the variable £. 

The two phases 9 and A must match to within an additive multiple of 27r, so 

2*o In (j^) = A(0 + 2vm, (8.9) 
where n is an integer. Using Eq. ([8.4)1 . we rewrite the relation in the form 

B 3 + — 2 = ^Me 2 ™/*o exp[A(£)/ So ]- (8.10) 
ma z m 

The values for n correspond to different 3-body bound states. Once the function 
A(£) is known, we are able to calculate the entire spectrum by solving Eq. ()8.10|) for 
various n. A(£) is what we have referred to as Efimov's universal function. 

The variables H and £ in our equations can be treated like a "radius" and an 
"angle" respectively. Suppose we plot a 3-body bound-state energy on a graph, 
choosing the x-axis to be 1/a and the y-axis to be — \J mB 3 /h 2 (See Fig. 18. 1|) . Then 
H is the distance of a line from that bound-state point to the origin, and £ is the angle 
that line makes with the x-axis. Figure 18.21 shows two Efimov states as a function of 
the 2-body binding energy which is proportional to 1/ma 2 to leading order. In this 
figure, we have chosen to label the axes differently so that the behavior of both states 
can be better seen. 

Using these parameters, we can demonstrate that the solutions to Efimov's equa- 
tion have a discrete scaling symmetry. If there exists a bound state with binding 
energy B 3 for the parameters a and A*, then there will also be a bound state with 
binding energy X 2 B 3 for the parameters A _1 a and A* provided that A is of the form 
exp [n'7i/s ] for some integer n' . This is equivalent to saying that for any given value 
of £, the values of H for any successive bound states differ by a multiplicative factor 
of exp(7r/so) — 22.7. Since we are considering only the case a > with three- 
body binding energies that satisfy B 3 > B 2 , the variable £ is restricted to the range 
— 7r/2 < £ < — tt/4. This periodic behavior can be easily seen by examining the value 
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of the phase in f d0 along a path of constant £. Figure 18.31 shows the value of 9 for 
bound states with B 2 /B 3 = 0.5 as H increases. Successive Efimov states must have 
the same phase which is periodic in \n(H). 

It should be noted that Eq. (j8.10|) is an approximation based upon the limit 
ro/a = 0, and as such is exact only for zero-range theories. Any calculated energies 
can be expected to have errors of O(r /a). Some first order corrections due to a 
non-zero effective range have been calculated by Efimov I3P] . 

8.1.1 Applications 

Efimov's general framework may be applied whenever there is a resonant two-body 
interaction. Several instances where it is applicable: 

• The two-nucleon system has a shallow bound state, the deuteron, in the spin- 
triplet channel and a large negative scattering length in the spin-singlet channel. 
This led Efimov to suggest that the few-nucleon system could be described using 
zero-range potentials with the effective range treated as a perturbation. The 
three-nucleon system includes two bound states: the triton (a pnn bound state) 
and 3 He (a ppn bound state). The equations are more complicated since the 
nucleons possess spin and isospin, and the ppn state is further complicated by 
the Coulomb interaction between the two protons. A leading order analysis 
with the Coulomb interaction neglected has been carried out by Efimov [TBI I3T] 
and revisited in the EFT framework by Bedaque, Hammer, and van Kolck |12j . 
An analysis at NLO in the effective range was carried out by Efimov j2H] and 
repeated by Hammer and Mehen using EFT [32J. 

• 4 He atoms have a large two-body scattering length and a shallow two-body 
bound state. The ground state of the 4 He trimer has been observed, but its bind- 
ing energy has not been measured. The Schrodinger equation for the three-body 
bound states has been solved accurately for potential models of the interaction 
between 4 He atoms. In addition to the ground state trimer, there is an excited 
state this is shallower by a factor of about 50 to 70. It has been found that 
Efimov's function can be used to predict the binding energy of one of the two 
3-body bound states using the other as input [121 ESI EH] • 

• For alkali atoms, the atom-atom scattering length can be made large by tuning 
and external magnetic field to a Feshbach resonance [221 122] ■ One complication 
in this case is that the alkali atoms have many deeply-bound two-body states, so 
an Efimov state can decay into an atom and a deep two-body bound state. As 
a consequence, the Efimov states are resonances with a binding energy B% and 
a width T 3 . The complex energies — (B 3 + ir 3 ) still satisfy Efimov's equation, 
but it requires the analytic continuation of the function A(£) to complex values 
of the angle £ [23]. 
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Level 


Energy 


Ratio 


+3 


1.366187266197138e+08 


515.0350013^ 


545557 


+2 


2.652610526516549e+05 


515.0350013^ 


550982 


+1 


5.150350013849171e+02 


515.0350013^ 


549171 





1.000000000000000 


515.0350013^ 


549775 


-1 


1.941615613134847e-03 


515.0350013^ 


551568 


-2 


3.769871 189167697e-06 


515.0350013^ 


546162 


-3 


7.319640760400369e-09 







Table 8.1: Efimov state energies for B2 = 0. The value 1.0 is fixed as a starting point. 
Notice that the ratio of energies for adjacent states equals exp(27r/so) to 11 digits. 



8.2 Computation of Efimov's Function 

Although Efimov explicitly considered only resonant two-body interactions, his con- 
clusions result from applying a boundary condition on the three-body wavefunction at 
short distances. This encompasses any short-range interactions, including the three- 
body contact interaction used in our work. The g% potential in our equations only 
acts when all three particles are very close together. Using our previous notation, this 
would correspond to the region R ~ r , which is inside the radius where the bound- 
ary condition matching occurs. Its influence is combined with that of the two-body 
interaction, and the total effect is seen only through the parameter A*. 



8.2.1 B 2 = 

As a verification of this argument, let us consider the case B2 = 0. According to 
Eq. ()8.10|) . the bound-state energies are given by 

B 3 = ^e 2 ™/ s « exp[A(-7r/2)/ So ], (8.11) 
m 

showing that the ratio of energies for adjacent states is exp(27r/so). To 13 digits, 
this ratio is 515.0350013848. In Table IHTT| we have computed several of the energies 
using B 3 = 1.0 as a reference point. The table also shows the ratios between adjacent 
states, all of which equal the predicted value to 11 digits. 

We can prove that the ratio of adjacent binding energies is exp(27r/s ) using 
Eq. (|7.1|) . Two adjacent bound states, B 3 and B' 3 , must have phases that differ by n. 
Therefore, 



s In (A//bT) = s In (a/V^T) + 7r 
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=> ^3=6^. (8.12) 
8.2.2 From to A 

Before we can compute Efimov's function, we must find a way to relate the parameter 
A* to the parameters in our equations. We will do this by considering the behavior 
of the phase 9 from the perspective of both the high- and low-momentum equations. 

First consider the high-momentum perspective. We have seen in Sec. I7.1l that it is 
possible to view 6 as a function of 8 , but it can just as easily be viewed as a function 
of G 3 . For a given value of G 3 , we write the phase generated by this coupling as 
h (G 3 ). 

From the low-momentum perspective, the phase is determined only by the param- 
eters B 2 , B 3 , and A. We also know explicitly the cutoff dependence of this phase: 



Oi(B 2 , B 3 , A) = s In \ k/\J B 3 j + 8(B 2 /B 3 ). (8.13) 
Because these phases must match to within a multiple of ir, 

6 h (G 3 ) = e t (B 2 , B 3 , A) + nix = s In (k/Jb~ 3 \ + 6(B 2 /B 3 ) + mr. (8.14) 



The purpose of G 3 is to make some value of B 3 cutoff-independent. As the cutoff 
is changed, G 3 is adjusted to ensure that B 3 remains the same. Removing A from 
Eq. ()8.14|) can only be done if 9h{G 3 ) contains A dependence of the form soln(A) 
implicitly through G 3 . We use this fact to define a new parameter A* via the equation 

e h {G 3 ) = s Q \n{K/K). (8.15) 

Eq. (J8.14)) can then be written as 

s In (A/A,) = s In (V/bT) + #(£ 2 /5 3 ) + nvr, (8.16) 

which implies 

B 3 = A2 e 2 ™ /so e 29 ~ /so . (8.17) 
The ratio B 2 /B 3 can be written in terms of £: 

B 2 /B 3 = l/tan 2 (-0. (8.18) 

Since 9 is a function of B 2 /B 3 , it is also a function of £, and from now on we will 
write #(£) to emphasize this. Multiplying both sides of Eq. ()8.17|) by (1 + B 2 /B 3 ) 
yields 



B 3 + Bo = A 2 e 2?m / S0 



exp 



0(O-21n(sin(-O) • (8.19) 
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We see that we can match Efimov's relation if we choose 

A* = Ae- e/so , (8.20) 
A(0 = ^(0-21n(sin(-0). (8.21) 

We should note that this choice for A* is proportional to the quantity A* defined in 
Ref. and will result in a function A(£) that differs from the one found in [22] by 
an inconsequential additive constant. 

8.2.3 Numerical Values 

Now that we have relations ()8.20|) and (|8.21|) . we can numerically compute values for 
A(£). One way to do this is to follow a state with constant A* as B 2 is changed. For 
a fixed cutoff, constant A* implies a constant 9. The two Efimov states in Fig. 18.21 are 
for constant phase, so they may be used to calculate A(£). The results are shown in 
Fig. 18.41 where we have plotted A(£) — A(— 7r/2) for both states. The constant shift 
in the function does not matter. 

A second alternative is to keep B3 constant and follow any changes in the phase 
as B2 changes. The function computed from this approach is shown in Fig. 18.51 for 
the case of B3 = 1.0. This method makes it easier to compute A(£) for any given £ 
since we do not need to search for the value of the three-body binding energy. 

Since the energy and phase values are accurate to about 12 digits, we expect 
similar accuracy in our calculation of A(£). Figure shows the difference between 
the Efimov function values generated from the Efimov state in Fig. l8.2l (-Bo = 0.0, B 3 = 
1.0) and the constant energy B 3 = 1.0 state. The errors support the statement 
that our calculation of A(£) is accurate to almost 12 digits, which is much higher 
than previous calculations [25 . Table I5~2l compares those previous values, labelled 
A(£)bhk, to the ones obtained here. The definition of these functions differ by an 
additive constant: A(£) = A(£)bhk + C- We have fixed the constant C by demanding 
that A(£) = A(£)bhk when £ = —1.502. The discrepancies are less than 0.02 with 
the exception of the final point at £ = —0.787 where the discrepancy increases to 
about 0.06. 



8.3 Number of Efimov States 

For any given value of B2, there will be a certain number of Efimov states that 
lie within the range of validity for this approximation. Efimov's estimate for the 
number of states N is related to the two-body scattering length a and the range of 
the interaction r by the simple formula 

N = - In (a/r ) . (8.22) 

7T 

This is easily derived from Eq. (|8.19J) . Assume a non-zero value for B2, which is 
related to the scattering length by B 2 oc a~ 2 . The deepest bound-state energy B 3 
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can only be of order A 2 before the conditions of our approximation are violated. This 
state will correspond to n = N, and we will write the energy as Bf ] ~ A 2 . The 
shallowest bound-state energy, corresponding to n — 0, must be greater than or equal 
to B 2 and should be of the same order. We write this as B$ ~ B 2 . Therefore, we 
have the relations 

B^ N) + B 2 = A 2 e 2 ^'exp(A(^)/ So ), (8.23) 
Bi 0) + B 2 = A 2 exp(A(£ )/s ), (8.24) 



where ^ = — arctan (yB^ jB 2 \ and£o = — arctan yy Bf^ /B 2 j. Since these states 

are part of the same spectrum, the value of A* in both equations must be identical. 

Because we are assuming that .83 ~ B 2 , we shall simplify Eq. ()8.24|) by making 
the substitution B^ 0) + B 2 ~ 2B 2 . Upon taking the ratio of (EHZHJ) and (JEHU), we 
obtain 

~ e 2 ^ exp [(A(6r) - A(£ )) /a ] ■ (8.25) 
The exponential will be 0(1), so we may write 

V ~ — In f — 77- I (8.20) 
2tt 2Bo ' 



The 2B 2 term in the logarithm can be approximated by a -2 . The cutoff A is related 
to the interaction range by A ~ rjT 1 , which allows us to estimate ~ Tq 2 . This 
makes our estimate of the number of states 



tf~^.ln [ a - =^ln(a/r ) (8.27) 
2vr V r o/ TT 

which matches the estimate given by Efimov. 

The source of this estimate can also be seen graphically. Recall Figs. 17.61 and 17.71 
where it is shown that successive bound states match onto the cosine behavior at 
successive peaks. The x-axis in these figures is labelled by x = s ln(p/A) so that 
the cosine behavior of fi as it approaches fdo can easily be seen. Each bound state 
joins the cosine curve at roughly x = soln^/A). If r]^ represents the shallowest 
state, then it joins the peak at x ~ s ln^f^/A). From Fig. 17.81 we see that the 
function fho decays exponentially when p ~ A. The last cosine peak is around xn — 
s ln(A/A) = 0. The number of peaks between xq and xn then gives the number of 
Efimov states that can exist: 

N = -(x N - x ) = -- hrf/A). (8.28) 

7T 7T 

Using A ~ Tq 1 and 773 ~ a -1 leads to the the same estimate N ~ (s /tt) ln(a/r ). 
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-6 


943 


-1 


214 


-7.15708283017 


-7 


151 


-1 


296 


-7.46636472602 


-7 


461 


-1 


347 


-7.63469982709 


-7 


632 


-1 


408 


-7.81593232651 


-7 


814 


-1 


443 


-7.91119231675 


-7 


910 


-1 


482 


-8.01059689781 


-8 


009 


-1 


502 


-8.05900000000 


-8 


059 



Table 8.2: Comparison of Efimov function A(£) to previously computed values. The 
values A(£)bhk are those computed by Braaten, Hammer, and Kusunoki j2Hj. The 
additive constant C was chosen so that the two functions agree at £ = —1.502. 
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Figure 8.1: Relation between bound-state energy and the variables H and £. 
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Figure 8.2: Binding energies of Efimov states as functions of the two-body binding 
energy. The dotted line is the scattering threshold. 
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Figure 8.3: Phase 9 for three-body bound states as a function of H with B 2 /B 3 = 
0.5. The periodic behavior as H increases is a direct result of the discrete scaling 
symmetry in Efimov's equation. For a constant ratio B 2 /B 3 , the phase is linear in 
ln(ff). 
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-0.5 -0.45 -0.4 -0.35 -0.3 -0.25 



Figure 8.4: The Efimov function A(£) as a function of £ computed using three-body 
bound states with constant phase. The difference A(£) — A(— 7r/2) is plotted since a 
constant shift in the function is inconsequential. 
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0.25 



Figure 8.5: The Efimov function A(£) as a function of £ computed using a three-body 
bound state with constant binding energy B 3 = 1. 
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Figure 8.6: Difference A err between the Efimov functions calculated using constant 
phase and constant energy as a function of £. The constant phase Efimov state with 
B 3 = 1.0 at B 2 = 0.0 is compared to the constant energy state B 3 = 1.0. The 
difference in the values suggests that our calculation is accurate to almost 12 digits. 



118 



CHAPTER 9 



CONCLUSION 



We have studied the three-body quantum mechanical bound state in the presence 
of short-range two- and three-body interactions. These short-range interactions allow 
us to study low-energy systems and compute results that are insensitive to the details 
of the inter-particle potentials. The equations describing the two- and three-body 
bound-state equations were derived using a regulated form of the delta function that 
removes divergences and lends itself to inclusion of higher order operators. 

To study the cutoff dependence of these equations perturbatively, we developed a 
method for uniformly expanding any function of three or four variables that possesses 
two widely-separated energy scales. This method was then applied to expand the 
three-body bound-state equation. Two sets of integral equations were derived: one 
for the leading order behavior, and one for the first order corrections. The solutions to 
these equations provide results that are accurate to 0(1) and 0{yfB^j 'A) respectively. 

We were able to derive several analytic results from the leading order equations, 
including a relation that allowed us to compare the parameters in our bound-state 
equation to those for Efimov's universal function. We proved the periodic behavior of 
the coupling 5q and used this to argue that G3 is periodic. The cutoff independence 
of all bound-state energies to leading order was also shown. 

Several techniques were shown that helped to improve the efficiency of our numer- 
ical computations. New limits on the integral equations were found, part of which 
involved substituting the analytic form for the asymptotic behavior in the low- and 
high-momentum regions. A method for achieving exponential convergence in the 
discretized integrals was given. 

The leading order equations were then solved numerically to high accuracy. We 
saw the limit-cycle behavior of G 3 , and displayed the divergent bound-state energies 
that make this coupling necessary. Solutions for the function were plotted to verify 
some of the assumptions we made in deriving the integral equations. We showed that 
our equation satisfied the conditions in Efimov's work, and proceeded to compute his 
universal function for three-body bound states to about 12 digits of accuracy. 
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There are a few areas available for possible future work. The integral equations 
for the leading order corrections have been derived, but they need to be solved nu- 
merically in order to study the leading cutoff dependence. The expansion method 
could also be taken further, and equations for 0(B 3 /A 2 ) terms could be derived. It 
may also be worthwhile seeing if the expansion method and numerical techniques 
used here can be of use in other non-related problems. 
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APPENDIX A 



ANALYTIC FORM OF MID-MOMENTUM 

FUNCTION 



This appendix derives the analytic form of fdi, 1 which satisfies the integral equa- 
tion 

, , v V [°° dq ( p 2 + q 2 + pq\ 

4n z D d0 {p) Jo q \p z + q z - pq j 

V f D dl {r] 2 ,p,A)\ po dq (p 2 + q 2 +pq\ 



47T 2 D d0 (p) V Aio(p) /Jo g \p 2 + q 2 -pq 
p ( D dl (r] 2 ,p,A)\ f°° dq f p 2 + q 2 + pq 



f*ln ?L±L±H (A.l) 



47r 2 73 d0 (p) V Aio(p) /Jo o \p 2 + q 2 -pq / 
with the -Dd functions reproduced for convenience: 

A»(P) = (A.2) 

l07T 



Ddi(V2,P,A) = VlV2 V3{h 2 + 8h 2 -16)p 
Ddo(p) V3p 8V2^A 

The key to solving this equation is understanding the solutions to 



(A.3) 



POO 

v{n,s,p)= I dq K{p,q)v{n,s,q) + p s [ln(p)] n , (A.4) 

JO 

where 

tf fe8 )= * J„(4±4±«). (A.5) 
V3 7r g \p 2 + q 2 -pqj 



1 This derivation is an expanded version of one provided by Wilson |15j . 
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The solution for n = is 



u(0,s,p) 



p 



i-A>(*)' 

8 sin (tt s/6) 
cos (its/ 2) 



(A.6) 
(A.7) 



except when s is an odd integer. For the special case of s = ±zsq, the function (3o(s) 
is one, and the solutions are then proportional to p s ln(p) [To] . 

Solutions for non-zero values of n can be obtained by differentiation. For example, 
taking the derivative d/ds of Eq. (jA.4|) for n = yields 



dz/(0, s,p) 
ds 



dq 



K(p, q) 



d z/(0, s, q) 
ds 



+ p s [Hp)} , 



(Ai 



which is the same equation satisfied by u(l, s,p). Thus, 

du(0,s,p) p s \n(p) 
^(1, = ; = tt-t^ + 



P 



d(3 Q (s) 



ds 



Po(s) 



Po(s)Y ds 



(A.9) 



The same procedure can be followed to obtain solutions for higher values of n. 

To make the connection to Eq. (|A.4|) clearer, the equation for fdi is rewritten as 



fdiip) 



+ 



+ 



dqK(p,q) f dl (q) 

V2t/2 \/3 (hj + 8h 2 - 16) p\ /■<» 



, \/3p 8V^vrA 

^(^ + 8/^2- 16) p\ 



, v^P 8\/2^A 
We will also make use of the fact that fdo satisfies the equation 



dq K{p, q) f d o(q) 
dqK{p,q)f dl {q). (A.IO) 



fdoip) 



dqK(p,q) f d0 (q), 



which has the solution 



f d0 (p) = A cos (s In (p/A) + 



(A.n; 



(A.12) 



A solution for fdi shall be constructed in a series of steps, ensuring self-consistency 
at each step. 

Consider Eq. (jA.10|) with only the first term on the right-hand side included. This 
equation is the same as the one for f d0 . Thus, fdi contains a term proportional to 
fdo- However, it must be of order rj/A, so we shall write this term as 



A d io -T- cos (s In (p/A) + dl0 ) 



(A.13) 
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where A d i is an arbitrary coefficient of 0(1), and 9 dw is an unknown phase. 
Using Eq. (jA.ll|) . the second integral in (jA.lO)) can be replaced by 



'V2V2 [hi + 8h 2 - 16) 



P 



fm{p). (A. 14) 



V V3p 8\/2^A 
If we added only the first of the two pieces, the integral equation for f d \ would become 

fdlip) = r dq K(p, q) f dl (q) + ^ f d0 (p). (A.15) 
Jo V3p 

We extend this equation into the complex plane by writing it as 

Ul{p) = f°° dq K(p, q) f dl (q) + ln ^ A )+ e ). (A. 16) 

Jo V3p 

The inhomogeneous term is proportional to p~ 1+1S0 ; so we may use Eq. ()A.4j) to obtain 
the solution 

™ = 1 - g „ ( - 1 + ,,) (S) ■ 4e "" , ~' < A ' 17 > 

The solution we desire is actually the real part of the previous solution. We shall 
write it as 

A u ^cos(s o hj(p/A) + n ). (A. 18) 

P 

The constants An and Q\\ are related to A and 9 by the relation 

(1 - /? (-l + is )) A n = J^-Ae ie . (A.19) 



If we had added the second piece of Eq. (IA.14J) instead of the first, our equation 
would have been 

U(p) = r d, K(p, ,) Mq ) - V3(W-16)P 
Jo 8v2vrA 

The steps for solving the equation for this case are nearly identical to the previous 
one. We extend the equation into the complex plane, noting that the inhomogeneous 
term is proportional to p 1+ls °. This solution is 

U(P) = ^ c ( ^ + 8fr 2 ~16M Aei(so Hp/A)+9) 

and we choose to write the real part as 

A hl | cos (so In (p/A) + B hx ) . (A.22) 
123 



The constants A h i and 6 hl are determined by 

j- ^(^ + 8^-16) 



;i-A)(l+iso)) A fcl e"« 



8V2^ 



Ae 



(A.23) 



At this point, we have three terms that contribute to fdi- Even though they were 
derived individually, they may simply be added together. The reason is that each term 
is "self-contained" in the sense that integrating it results only in a term proportional 
to itself. There are no cross-terms generated in the equation. The solution, as it now 
stands, is 

/di(7/2,p,A) = A a — cos(s hi(p/A) + a ) 

+ A hl jCos(s ln(p/A) + 9 h i) 

+ A dl0 j cos (s In (p/A) + d io) • (A.24) 

The only thing remaining that affects the solution is the third integral from 
Eq. (jA.lOJl . Using the version of fdi from (|A.24|) . its contribution is 



' V2 V2 _ y / 3(^ + 8/i 2 -16)p x 
, V3p 8^ A 

fV2r] 2 VS(h 2 2 + 8h 2 -lQ)p 



dqK{p,q)f dl (q) 



V2 



V 



8V27T A 



A d \v cos (s In (p/ A) + dlo ) 



+ 



+ 



cos (s In (p/A) + 9a) - ^ /^(p) 
p V3p 

P / i / l \ \ n \ ^>{h 2 + 8h 2 - 16) p , 
Awy cos (s In (p/A) + 9 hl ) + ~ r— ^ / d0 (p) 



VS(h 2 + 8h 2 - 16) 



8V2^A 



An cos (s In (p/A) + 



+ £^ Ai cos (s In (p/A) + fc i) 
(/il + 8h 2 - 16) ?7 2 



4v / ttA 



A cos(s ln (p/A) + 0) . 



(A.25) 



Note that we have dropped all terms that are 0(rf / A 2 ) or higher. Although this 
contribution can easily be absorbed into the A^io term, we must also consider the 
terms it may generate by acting as an inhomogeneous term in our integral equation. 

To keep the equations simple, we will combine these three terms into one. Since 
we will be extending them into the complex plane, let us define 



Hi p^ S0 ln (p/A)+e dl ) 
A 



VS {h 2 2 + 8h 2 - 16)_?72 A e i isQ i n{p/ A )+ e n ) 



8V2^A 
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+ V2j]2 A e i (so i„ (p /A)+e hl ) 



V3A 



{hi + 8h 2 - 16) 7/2 



^ e i(s ln(p/A)+^ 



40FA 

We must now solve the equation 

fdlip) = f°° dq K(p, q) f dl {q) + ^ A dl '^'^l 
jo A 

The solution takes the form 

d 



fdiip) 



ds 



A)00 



In 



S=1S 



^ !^ 4 Jt pi(^o ln(p/A)+O dl ) 

A/ A 



(A.26) 



(A.27) 



(A.28) 



The real part will be written as 



A dU ^ In (p/A) cos (s In (p/A) + 0*) 



subject to the relation 



^irfiie — ^i<iie 



(A.29) 



(A.30) 



-WO 



Equation ()A.29|) adds a new term to the solution. The only remaining considera- 
tion is whether or not it affects the results of Eq. (jA.25J) since the value of fdi used in 
that equation did not contain this new term. Fortunately, the new term would only 
contribute at 0{rf / A 2 ) so the result remains unaffected. 

We now have the self-consistent, analytic solution shown previously in Chapter |S| 

fdi(V2,p,A) = 4i— cos (s In (p/A) + 0u ) 
p 

+ A hl |cos(s ln(p/A) + 9 hl ) 
+ A diojr cos (s In (p/A) + 9 dw ) 

+ Adnj In (p/A) cos (s In (p/A) + 6 dn ) . (A.31) 
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APPENDIX B 



PROOF OF EXPONENTIAL CONVERGENCE 



The proof given here is based on notes provided by Kenneth Wilson. It is meant to 
give the reader a justification as to why the uniform spacing of integration points on 
a logarithmic scale gives exponential convergence, but it will not be entirely rigorous. 

We start by making the following assumptions: 

• f(z) is analytic in the region — oo < Ke(z) < oo and —r < lm(z) < r for some 
value of r. 

• f(x) is real for any real value x. 

• The integral /f^ f(x) dx is finite. 

The first two assumptions imply (via the reflection principle) that f(z*) = f(z)* in 
this region With this in mind, consider 



dx 



f(x + ir) 



1 — exp(— 27ri(x + ir) /a) 



+ 



fix 



ir 



exp(27ri(x — ir)/a) 



(b.i; 



This must be real since the second term is the complex conjugate of the first. There- 
fore, 



dx 



f(x + ir) 



< 



CO 

-oo 
oo 



1 — exp(— 27ri(x + ir)/a) 
f(x + ir) 



+ 



fix 



17" 



dx 
dx 



exp(27ri(x — ir)/a) 
f(x - ir) 



+ 



1 — exp(— 27ri(x + ir)/a) 1 — exp(27ri(a; — ir)/a) 
f(x + ir) f(x — ir) 



+ 



< e -2wr/a 

g— 2nr/a 



1 — exp(— 27ri(x + ir)/a) ' 1 
f(x + ir) 



exp(27ri(x — ir)/a) 



dx 



< 



-2-irr/a 



exp(— 2-Kr/a) — exp(— 27ria;) 
dx \f(x + ir) \ + \ f(x — ir) \ . 



+ 



f(x- ir) 



exp(— 2irr/a) — exp(— 27ria;) 

(B.2 
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However, if f(z) is such that it tends to zero as \z\ — > oo, then from the analyticity 
of f(z) we know that 



dx f(x) 



dx f(x + ir) 



dx f{x — ir). 



(B.3) 



Now 



/OO rOO rOO 

<ix |/(x + ir)| + |/(x — ir)| < / dx Re[f(x + ir)] + / dxlm[f(x + ir)] 
-oo J — oo J— oo 

/oo /*oo 
<ix Re[/(x — ir)] + / <ix Im[/(x — ir)] 
-oo J —oo 

/OO /"OO 
dxRe[/(x + ir)] + / dxRe[/(x — ir)] 
-oo J —oo 

dxf(x). (B.4) 



oo 
oo 



Combining this with Eq. (jB.2|) yields 

f(x + ir) /(x — ir) 



dx 



^ g— 27ri(x-|-ir)/a J g2-7ri(x— ir)/a 



< 



g— 2-nr/a 



dxf(x). (B.5) 



If exp(— 2'Kr/a) is small, then the above term is small provided that the integral is 
not unusually large. We use this to write 



dx f(x) 



dx f(x) — 



dx 



f(x + ir) 



+ 



f(x - ir) 



J g— 27ri(a;+ir)/a ^ g27ri(x— ir)/a 



(B.6) 



The second and third terms have poles on the real axis at na where n — 0, ±1, ±2, 
etc. This means we can move their contour integrations from x + ir and x — ir to 
x + ie and x — ie respectively (where e is some small quantity) without crossing any 
poles. Therefore, Eq. (jB.6|) becomes 



dx f(x) 



dxf(x) — j dx 



f{x + ie) 



+ 



/(x - ie) 



~^ g— 27ri(x-|-ie)/a ^ g2-7ri(a;— ie)/a 



(B.7) 



Now consider a very small loop around the entire real axis. The Residue Theorem 
implies that 



°° dx fj£ - ie ) 
-oo 1 - e 27rl ( :r - le )/ a 



°° /(x + ie) . " a/(na) 

ax Q ., ,. w = 27ri > 

1 p 27ri(x+ie)/a 

1 n=— oo 



oo 



2vri 



(B.8) 



where — a f (an) /2ni is the residue from the pole at x = na. Combining Eq. (jB.7|) and 
Eq. (|B.8|) gives us 



dx f(x) 



dx /(x+ie) 



1 



1 



^ g— 2-7ri(x+ie)/a ^ g27ri(a;-)-ie)/a 



a f(an). 

i=—oo 

(B.9) 
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Here we have replaced / dx f(x) with / dx f(x + ie) which makes a small error of 0(e). 
The term in brackets is equal to zero, leaving us with 



/oo °° 
dxf{x) ~ o/W. (B.10) 
-oo 



The error in this approximation is equal to the magnitude of Eq. (jB.5|) . which goes 
like exp(— 2irr/a). The parameter a represents the spacing used in discretizing the 
integral. As it approaches zero, the error decays exponentially. 
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APPENDIX C 



ERROR ANALYSIS 



This appendix provides a simple error analysis for the integration technique we 
use. It is meant to justify the claim of exponential convergence. If the answers truly 
converge exponentially, then we expect the approximation using n integration points 
to behave like 

A n = A + ae~ pn , (C.l) 

where A represents the exact value. By taking some large value n = N and using A^ 
as the value for A, we should find that 

e „ = io g J^^U-^»M n + !^, (c.2) 

\ A N ) A N A N 

where e n is the relative error in A n . Plotting e n as a function of n should show a 
straight line. At some value of n, other sources of error, such as machine error, will 
dominate and e n will become roughly constant. The value of e n in this flat region 
tells us the number of digits of accuracy in our result. 

The integral equation for f i0 was used to calculate the phase using different num- 
bers of integration points. This was done for three sets of B2 and B3 values, and 
the results for the relative errors are shown in Fig. 10. 1 1 The initial linear behavior 
indicates exponential convergence, which is then dominated by other errors resulting 
in a plateau. The plot indicates an accuracy of about 13 digits. 
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Figure C.l: Relative error e n in the phase calculation for different values of B 2 and B 3 
as a function of the number of integration points n. Three different sets of energies 
were used. Notice the initial linear behavior indicating exponential convergence, and 
the plateau indicating the limit of accuracy (approximately 13 digits). 
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